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FOREWORD 


This  document  describes  the  work  performed  by  the  Philco-Ford  Corporation  under 
the  Advanced  Penetration  Problems  project  (APP)  Contract  F0470 1- 7 L-C-0035  in 
the  period  1971-1972.  In  the  present  period,  emphasis  was  placed  on  the  study 
of  boundary- layer  turbulence  at  supersonic  and  hypersonic  Mach  numbers,  on  the 
effect  of  altitude  changes  on  the  structure  of  the  near  wake  behind  a  cone, 
and  the  statistical  features  of  far  wakes  and  ionized  Jets.  The  principal  in¬ 
vestigator  of  the  project  Is  Dr.  Anthony  Deinetrindes  of  the  Fluid  Mechanics 
Section,  Research  Laboratory,  at  the  Newport  Beach  facility  (Telephone  640-1500, 
Extension  1841).  Project  Scientists  are  Dr.  Demetriades  and  Dr.  Arnold  Laderman, 
also  of  the  Fluid  Mechanics  Section  (Telephone  640-1500,  Extension  1012).  The 
work  has  had  the  able  support  of  Mr.  L.  Von  Seggern  and  Mr.  E.  L.  Doughman  at 
the  Philco-Ford  Supersonic  Wind-Tunnel. 

The  present  effort  was  sponsored  jointly  by  SAM, SO,  United  States  Air  Force  and  the 
U.  9.  Army  Advanced  Ballistic  Missile  Defense  Agency  (ABMDA) .  In  addition  to  the 
Philco-Ford  tunnel, facilities  utilized  for  tests  included  the  21"  Hypersonic  Wind- 
Tunnel  at  NASA's  Jet  Propulsion  Laboratory  and  TunneL  B  at  the  Von  Harman  Facility 
at  AEDC,  Tullahoma,  Tennessee.  The  courtesy  extended  by  the  General  Electric  Co. 
in  making  available  the  porous  cone  model  for  the  latter  tests  is  gratefully  ack¬ 
nowledged.  Also,  the  cooperation  is  acknowledged  of  the  Von  Karman  Facility  per¬ 
sonnel  for  designing,  fabricating  and  installing  the  hot-wire  shield  and  actuator 
device  utilized  at  Tunnel  B. 

This  report  has  been  reviewed  and  approved. by 

Leonard  T.  Atkins 
Contracting  Officer 

Space  and  Missile  Systems  Organization 
Mr  Force  Systems  Command 
Los  Angeles,  California  90009 
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ABSTRACT 


Thu  study  of  turbulence  on  a  flat  plate  at  Mach  9.4  has  been  completed  with  the 
discovery  that  the  pressure  fluctuations  predominate,  that  the  velocity  fluctua¬ 
tions  are  very  small  and  that  wall-cooling  has  a  strong  influence  in  decreasing 
the  temperature  fluctuations.  Tests  have  been  performed  at  AEDC/B  to  extend 
these  results  to  the  cone  geometry,  to  broaden  the  Mach-  and  Reynolds-numbers 
matrix  and  to  investigate  ablation  effects.  Although  still  incomplete,  analysis 
of  the  data  show  very  largo  temperature  fluctuations  due  to  the  adiabatic  wall, 
relatively  little  change  of  the  turbulence  level  past  the  transition  point,  and 
discernible  but  small  fluctuation  increase  for  small  mass  injection  rates.  Cor¬ 
relations  have  began  of  all  the  available  boundary- layer  turbulence  data  for  the 
purpose  of  scaling  the  fluctuation  distribution  to  arbitrary  flight  conditions. 
Thus  far  the  available  correlation  schemes  fit  the  velocity  fluctuation  data 
only  in  the  outer  half  of  the  layer.  Better  success  is  enjoyed  for  the  temper¬ 
ature  difference  between  the  wall  and  the  external  flow.  Experiments  were  also 
conducted  to  measure  the  transverse  velocity  fluctuation  component  and  the  tur¬ 
bulent  shear  stresses  with  specially-developed  hot-wire  probes;  these  data  are 
now  available  on  tape,  In  the  near-wake  investigation  the  completed  mean-flow 
measurements  were  analyzed  and  show  that  a  "plateau"  tends  to  form  as  the  Reynolds 
number  increases  beyond  BLT.  The  turbulent  viscosity  was  measured  and  was  found 
to  decrease  beyond  BLT  as  predicted  by  the  AVCO  theory.  Other  theoretical  pre¬ 
dictions,  such  as  the  radial  distribution  of  velocity  beyond  the  neck,  were  also 
verified.  Fluctuation  measurements  reve.aled  that  the  shoulder  expansion  nearly 
destroys  the  velocity  but  not  the  temperature  fluctuations,  and  that  the  wake 
is  "double-structured"  beyond  the  neck.  The  study  of  electron,  velocity  and 
temperature  space-time  correlations  in  the  far  wake  has  been  completed. 
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Vertical  displacement  between  the  bottom  of  the  surface  sense  probe 
and  the  probe  i 

D.C.  voltage  output  of  wide  band  data  (wake  test) 
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Non-dimensional  narrow-band  fluctuations  j 

Noise  equivalent  of  Z  '>j 

Function  of  Mach  number  1 

First  hot-wire  resistivity  coefficient  at  0°C  ’ 

i 

Function  of  Mach  number  J 

Ratio  of  specific  heats 

Boundary  layer  thickness  ,-j 

Sublayer  thickness 

Drag 

Recovery  factor;  also,  coordinate  'y/h  "] 


NOMENCLATURE  (Continued) 


0  Momentum  thickness;  also,  temperature  defect;  also,  cone  half-angle 

0 (i)  Angle  between  vertical  and  straight  line  connecting  probe  i  to  cone 

center 

A  Integral  turbulence  scale 

X  Non-dimensional  blowing  rate 

v  Kinematic  viscosity 

“T  Eddy  kinematic  viscosity 

E  '  Instantaneous  transverse  component  of  vorticity  fluctuation  normal¬ 

ized  with  respect  to  local  mean  vorticity 

rr ’  Instantaneous  sound  field  fluctuation  normalized  with  respect  to 

strength  of  local  mean  sound  field 

density  fluctuation  normalized  with  respect  to  local 

entropy  fluctuation  normalized  with  respect  to  local 

Turbulent  Prandtl  number 

T  (T'V 

t'  Instantaneous  longitudinal  component  of  vorticity  fluctuation  normal¬ 

ized  with  respect  to  local  mean  vorticity 

T  Wall  friction 

w 

0  Angle  between  hot-wire  element  and  mean  flow  velocity  vector;  also, 

circumferential  angle;  also,  slope  of  Nu-Re  curve 

i|i  Slope  of  recovery-factor  versus  Mach  number  curve 


0  Density 

p'  Instantaneous 

mean  density 


o'  Instantaneous 

mean  entropy 


xv  ii 


NOMENCLATURE  (Continued) 


SUBSCRIPTS 


e  Edge  condition 

o  Stagnation  conditions 

r  Recovery  conditiona 

rms  Root  mean  square  value 

w  Wall  condition 

x,y,z  In  the  x,  y,  z  directions 

«  Free-Btream  conditions 


SUPERSCRIPTS 

(  )  Denotes  time  averaged  mean  value 


(  )  ' 


RMS  fluctuation 


SECTION  I 


SUMMARY  OF  MAJOR  FINDINGS 


In  the  present  contract  period  a  major  portion  of  the  work  was  concentrated  on 
the  definition  of  hypersonic  boundary- layer  turbulence  properties.  In  actual 
flight  the  problem  requires  knowledge  of  a  thin  boundary- layer  (6  «  the  local 
body  section  diameter)  on  a  cone  with  external  Mach  number  Me  a  10  to  12,  Reynolds 
number  Rg  from  10^  to  10^,  heat  rates  such  that  0.05  <  Tw/T0  <  1  and  ablation 
rates  m/ceue  ranging  from  0  to  0.01,  Since  it  is  nearly  impossible  to  produce 
these  conditions  in  a  single  ground  test  suitable  for  turbulence  measurements, 
scaling  laws  are  required  which  can  be  formed  by  tests  over  a  wider  matrix  of  con¬ 
ditions  . 

Data  from  the  Mu  -  9.4,  Rg  -  37,000  tests  at  ' PL/HWT  (code-named  "TBL-A")  com¬ 
pleted  during  1971-1972  were  re-examined  and  lined  (see  Section  II),  and  a  sum¬ 
mary  Technical  Report  has  been  Issued  separately  (Reference  1).  Doubts  cast 
earlier  on  the  ability  of  this  test  to  represent  highly  "mature"  layers  have  been 
dispelled  because  the  vertical  gradient  5p/dy  has  since  been  found  to  be  independent 
of  the  wind  tunnel.  The  pressure  fluctuations  p'  have  been  found  to  be  the  dom¬ 
inant  turbulence  mode  at  high  Me,  reaching  18%  intensity,  with  the  density,  tem¬ 
perature  and  velocity  fluctuations  following  in  that  order.  Intetmittency  measure¬ 
ments  made  from  the  TBL-A  tape  library  also  showed  an  unexpectedly  smooth  ("un¬ 
wrinkled")  turbulent  front  implying,  for  example,  reduced  radar  scattering  from  a 
very  overdense  layer  as  opposed  to  what  would  be  predicted  from  low-Mach-number 
boundary- layer  experience. 

To  apply  the  flat-plate  lessons  learned  at  JPL  to  a  flight-realistic  cone  geometry, 
to  supply  further  variations  in  M„  and  Rg  and  to  introduce  ablation  simulation, 
tests  have  been  performed  at  AEDC/B  (code-named  "TBL-C")  in  the  experiments  dis¬ 
cussed  in  Section  III.  With  the  cooperation  of  AEDC  personnel  the  initially  ser¬ 
ious  problem  of  hot-wire  anemometer  survival  and  endurance  in  Tunnel  B  was  re¬ 
solved  successfully.  The  tests  were  then  performed  at  Mach  8  (Me  ~  7)  at  various 
combinations  of  Rg,  stations  x  from  the  cone  tip  and  mass  injection  rates.  The 
data  are  presently  under  reduction,  but  sufficient  preliminary  information  is 
available  to  draw  some  conclusions.  For  example,  it  appears  that  mass  injection 
as  high  as  l  s  m/peue  -  0.0015  does  not  affect  the  maximum  fluctuation  magnitude 
in  the  layer.  It  also  appears  that  the  relative  distribution  of  the  velocity-vs- 
temperature  fluctuations  is  not  changed  at  AEDC  over  what  had  been  found  at  JPL. 

With  such  information  now  becoming  available,  a  separate  effort  was  set  up  to  col¬ 
lect  and  correlate  the  turbulence  data  from  TBL-A,  TBL-C  and  other  sources  into  a 
coherent  picture  (Section  IV).  Increasing  Me  was  found  to  cause  a  decrease  of  the 
velocity  fluctuations  u'.  Tentatively,  n  correlation  of  the  form  (u'/uT) (p/pw)* 
vs  y/6  was  found  to  describe  all  data  fairly  well  for  y/6  >  0.5,  but  only  poorly 
closer  to  the  wall,  for  a  wide  range  of  Me,  q  and  Rg.  On  the  other  hand,  the*  tem¬ 
perature  fluctuations  T'  increase  rapidly  with  Me,  appear  to  be  independent  of  Rg, 
and  decrease  rapidly  as  the  wall  temperature  decreases  below  the  adiabatic-recovery 
level.  Fiir- to-good  success  lias  been  found  in  correlating  the  T1  data  with  Kist lev's 
original  suggestion  that  the  "driving  force"  for  T'  is  the  mean  temperature  drop 


Tw  -  Te  across  the  layer.  Finally,  Che  longitudinal  integral  scales  of  all 
variables  but  Che  pressure  appear  to  follow  a  geometrical  relation  A  —5/5  inde¬ 
pendent  of  Me,  Rg  and  heat  rate. 

To  measure  the  transverse  turbulent  components  v'  and  the  turbulent  shear  stresses 
u'v1  and  p 1 v  1  ,  a  test  code-named  TBL-B  was  performed  during  this  period  with  spec¬ 
ially-made  crossed  hot-wires  ("X  probes")*  To  reduce  time  and  costs  the  now- 
familiar  boundary- layer  at  JPL  was  again  used.  Section  V  describes  the  tests; 
this  measurement  was  the  first  attempt  ever  made  to  obtain  directly  the  ingredients 
of  the  "eddy  viscosity"  at  speeds  other  than  Me  ■  0.  By  also  obtaining  v'  it  will 
be  possible  to  resolve,  at  least  partly,  the  question  of  isotropy  which  is  usually 
assumed  to  hold  in  the  boundary- layer .  Provisions  to  reduce  these  data,  which  now 
exist  on  tape,  had  not  been  made  in  the  contract;  data  reduction  is  provided  for, 
however,  in  the  1972-1973  period. 

The  boundary- layer  turbulence  also  defines  initial  conditions  for  wake  formation. 
The  wake  development  subsequent  to  the  cone  shoulder  formed  the  second  major  ef¬ 
fort  area  under  this  contract.  Previously,  a  5-degree  half-angle  cone  had  been 
placed  in  the  Philco-Ford  Supersonic  Wind-Tunnel  (PF/SWT)  and  its  wake  subjected 
to  mean- flow  measurements  in  an  experiment  code-named  WEM.  In  the  present  period 
turbulence  data  were  acquired  behind  the  Bame  cone  as  described  in  Section  VI. 

The  tests  are  performed  in  a  background  of  urgent  discussion  in  the  defense  com¬ 
munity  about  the  peculiarities  of  a  wake  formed  when  the  vehicle  boundary- layer  is 
turbulent  ("below  BLT").  The  main  contention  is  that  the  shed  boundary- layer 
structure  decreases  the  wake  eddy  viscosity  by  its  small  Hcale  and  causes  a  delay 
in  wake  diffusion.  However  there  had  previously  been  no  conclusive  evidence  that 
this  contention  was  true.  During  this  period  continued  analysis  of  the  mean-flow 
data  in  (J  <  x  <  6D  (0  <  x  <  45(CjAp)  showed  that  an  increase  of  flight  Reynolds 
number  Rep)  towards  BLT  and  beyond  causes  an  increase  and  a  "flattening"  of  the 
wake  velocity.  Furthermore  the  eddy  viscosity,  when  extracted  from  the  data, 
showed  the  claimed  decrease  from  the  pre-BLT  to  the  post-BLT  case,  On  the  other 
hand,  the  seemingly  contradictory  assumption  of  a  Gaussian  radial  velocity  pro¬ 
file  beyond  the  neck  was  also  verified.  In  fact  these  profiles  are  numerically 
almost  identical  to  the  known  equilibrium  far-wake  profiles. 

Quantitative  turbulence  measurements,  covering  most  of  the  required  matrix,  were 
also  made  during  this  period.  The  technique  used  measures  both  the  fluctuation 
intensities  and  their  scales,  from  the  body  shoulder  to  40(CjjA)'«  .  It  was  found 
that  the  shoulder  expansion  indeed  decreases  the  velocity  fluctuations  u1,  but 
that  it  hardly  affects  T'.  This  provides  an  explanation  of  the  "turbulence  per¬ 
sistence"  seen  in  ballistic-range  photographs.  Indeed,  due  to  the  "Strong  Reynolds 
Analogy"  (i.e.,  kinetic-thermal,  energy  exchange)  the  near-wake  consists  of  a  core 
of  velocity  fluctuations  surrounded  by  an  annulus  of  temperature  fluctuations. 
Presently  the  shear  layer  and  the  boundary- layer  are  found  to  be  equally  respon¬ 
sible  for  creating  this  double  structure. 

In  addition  to  the  near  wake  studies,  two  projects  remaining  from  an  earlier  con¬ 
tract:  period  have  been  concluded.  The  space-time  correlations  ol'  an  equilibrium 
("far")  wake  have  been  mapped  and  are  now  available  for  numerical  insertion  into 
scattering  calculations.  Also,  the  functional  relations  have  been  obtained  be¬ 
tween  electron  and  gas-dynamical  fluctuations  in  a  turbulent  plasma.  Separate 
reports  on  these  two  projects  will  be  issued  soon, 


In  the  course  of  the  AEDC  tests  at  M  ■  8,  It  was  observed  that  for  a  consider* 
able  distance  preceding  transition  tne  boundary-layer  ahowed  Intense  sinusoidal 
instabilities.  These  waves  have  been  previously  observed  photographically,  but 
never  before  quantitatively  with  turbulence-measuring  probes.  The  wavelength 
found  is  approximately  1.76,  the  intensity  nearly  30'/.  (a  large  number  consider¬ 
ing  that  the  turbulence  levels  were  found  to  be  about  the  same)  and  the  spectral 
composition  almost  purely  temperature.  Since  similar  waves  photographed  else¬ 
where  show  the  same  wavelength,  it  has  been  concluded  that  this  is  a  real  effect 
which  is  important  to  the  transition  process. 
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SECTION  II 


TURBULENT  BOUNDARY  LAYER  MEASUREMENTS  AT  M  -  9.4 


2.1  SUMMARY  OF  PREVIOUS  WORK 

Under  the  APP  project,  Philco-Ford  has  been  investigating  in  detail  the  mean  and 
fluctuating  flow  in  a  highly  equilibrated  C Rq  “  36,000)  turbulent  boundary-layer 
over  a  cooled  wall  at  M®  ■  9.4.  The  nozzle  wall  of  the  Jet  Propulsion  Laboratory' 
Hypersonic  Wind-Tunnel  (JPL/HWT)  has  been  utilized  because  of  its  low  cost,  ease 
of  access  and  the  inherent  advantages  to  the  diagnostic  methods  of  the  large  thick 
neas  (4  inches)  of  its  boundary- layer .  This  project  code-named  "TBL-A"  is  the 
first  serious,  quantitative  attack  on  the  high-speed  boundary- layer  turbulence  pro 
blem  since  1957  and  the  first  ever  into  its  truly  hypersonic  regime, 

During  the  1971-1972  period  the  mean  flow  tests  on  this  flow  were  completed  and 
the  results  presented  in  Reference  2.  For  background  purposes  these  results  are 
abstracted  below  from  that  reference! 

"An  experimental  investigation  was  carried  out  to  determine  the  structure 
of  the  turbulent  boundary-layer  on  the  wall  of  the  Jet  Propulsion  Labora¬ 
tory  Hypersonic  Wind-Tunnel  at  a  free-stream  Mach  number  of  9.37.  Pro¬ 
files  of  flow  properties  were  obtained  from  pitot  pressure,  static  pres¬ 
sure,  and  total  temperature  surveys  made  through  the  4-inch  thick  boundary- 
layer  at  a  station  160  inches  from  the  nozzle  throat.  Tests  were  conducted 
primarily  for  the  following  conditions!  unit  free-stream  Reynolds  number 
of  127,000  per  inch,  corresponding  to  a  Reynolds  number  based  on  momentum 
thickness  of  36,800,  and  a  wall-ta-free-stream  total  temperature  ratio  of 
0.385.  A  cursory  examination  of  the  boundary  layer  was  also  made  at  a 
unit  stream  Reynolds  number  of  67,000.  The  static  pressure  measurements, 
although  indicating  negligible  streamwi.se  variation,  revealed  the  exist¬ 
ence  of  a  significant  pressure  gradient  normal  to  the  wall,  with  the  pres¬ 
sure  at  the  wall  approximately  45%  greater  than  its  free-stream  value, 

The  data  indicate  that  the  profile  of  total- temperature-ratio- versus- 
velocity-ratio  agrees  with  the  linear  Crocco  theory  in  the  sub-layer  region, 
but  in  the  outer  portion  of  the  boundary-layer,  closely  follows  the  quad¬ 
ratic  relation  which  has  been  observed  to  characterize  norszle  wall  measure¬ 
ments.  The  data  also  indicate  that  the  sub- layer  is  0.05  to  0.10  inches 
thick  which  corresponds  to  only  several  percent  of  the  total  boundary-layer 
thickness.  Correlation  of  the  velocity  profile  with  the  conventional  in¬ 
compressible  profile  shows  poor  agreement  in  the  sub-layer  and  wall-of  the- 
wake  regions,  with  the  latter  attributed  primarily  to  the  pressure  gradient 
normal  to  the  wall.  However,  the  experimental  skin  friction  coefficient 
agrees  with  the  value  predicted  by  the  Van  Driest  theory  within  several 
percent.  Finally,  the  total  temperature  profile  determined  from  a  hot¬ 
wire  anemometer  traverse  through  the  boundary- layer  is  in  good  agreement 
with  the  thermocouple  measurements." 

The  results  of  the  turbulence  study  also  completed  to  a  large  extent  in  the  1971- 
1972  period,  were  issued  under  a  separate  report  (Reference  1)  from  which  the 
abstract  is  produced  below: 


"An  experiment  has  been  performed  to  measure  the  distribution  of 
fluctuations  in  a  turbulent,  cooled-wall  boundary- layer  in  contin¬ 
uous  flow  at  Mach  9.4.  The  0.00001  inch  diameter  constant-current 
hot-wire  anemometer  was  employed.  The  hot-wire  signals  havu  been 
interpreted  both  with  and  without  aid  from  the  "no-sound"  assumption. 

Both  signal  interpretation  schemes  give  results  within  307.  of  each 
other.  The  fluctuation  intensity  in  the  boundary-layer  was  found  to 
consist  mainly  of  high-frequency  pressure  fluctuations,  whose  magni¬ 
tude  at  the  wall  and  beyond  the  layer  edge  agree  with  expectations. 

Static  temperature  fluctuations  are  lower  than  expected,  and  it  is 
suggested  that  this  is  caused  by  wall  cooling.  The  fluctuations  in 
the  longitudinal  velocity  component  are  generally  small  and  numeri¬ 
cally  differ  little  from  the  lower  Mach  number  results.  The  maturity 
of  the  boundary- layer  was  judged  high  because  of  the  measured  large 
turbulence  Reynolds  number  and  the  observed  inertial  subrange  In  the 
spectra.  In  the  latter  no  changes  were  noted  from  one  convection 
variable  to  another,  but  the  pressure  spectra  are  different  from  the 
velocity  and  temperature  spectra.  The  sublayer  and  boundary- layer 
thicknesses  were  determined  by  intermittency  measurements,  which  also 
showed  that  turbulence  bursts  occur  throughout  the  sublayer  and  oc¬ 
casionally  reach  the  wall.  The  standard  deviation  of  the  boundary- 
layer  front  is  much  smaller  than  observed  at  low  speeds." 

Additional  points  of  Importance  have  been  brought  to  light  during  the  present 
period;  these  are  described  below. 

2.2  FURTHER  TURBULENCE  CHARACTERISTICS  OF  THE  BOUNDARY- LAYER 

2.2.1  CONSISTENCY  CHECKS!  PRESSURE  GRADIENT  AND  SUBLAYER  THICKNESS 

Since  this  experiment  was  crucial  because  it  produced  the  long-desired  turbu¬ 
lence  information,  it  is  important  to  ensure  that  its  boundary- layer  was  not 
limited  in  its  applicability.  For  example,  we  had  observed  a  pressure  gradient 
dp/ay  such  that  the  wall  static  pressure  pw  ■  1.43  p*,  (see  Figure  15,  Reference 
2).  This  was  considered  unusual  and  was  at  first  met  with  some  skepticism; 
suggestions  were  advanced  Chat  it  was  peculiar  to  dynamic  non-equilibrium  ef¬ 
fects  in  wind-tunnel  nozzles.  Since  then,  however,  Owen  (Reference  3)  has  pre¬ 
sented  a  correlation  of  a  similar  effect  for  a  number  of  experiments  in  the 
range  2  <  <  40,  Significantly  his  correlation,  reproduced  here  in  Figure  L, 

includes  data  on  various  bodies  such  as  cones,  wedges  and  flat  plates  in  addi¬ 
tion  to  data  from  nozzle  walls.  The  Philco-Ford  point,  marked  "TBL-A"  ia  not 
far  from  Owen's  correlation.  The  conclusion  is  that,  within  experimental  ac¬ 
curacy,  the  pressure  gradient  found  by  Philco-Ford  at  JPL/HWT  is  a  true  effect 
not  limited  to  wind-tunnel  nozzle  waLls. 

Another  criterion  for  normalcy  is  the  sublayer  thickness  6 l .  Figure  2  shows 
another  finding  by  Owen  (Reference  3)  in  which  &l/6  is  correlated  with  the  para¬ 
meter  Me/Rj}^.  Again,  the  Philco-Ford  results  are  in  satisfactory  agreement  with 
this  correlation,  especially  since  the  definition  of  6 ^  and  6  are  subject  to  some 
error.  In  our  case  we  chose  6^  as  found  from  the  intermittency  data  (see  Section 
III,  Reference  1). 
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FIGURE  1.  EVIDENCE  OF  PRESSURE  GRADIENT  NORMAL  TO  THE  WALL  AT 
SUPERSONIC  AND  HYPERSONIC  SPEEDS . 


2.2.2  STATUS  OF  TBL-A  RESULTS:  NEW  FINDINGS 


Following  some  minor  refinements,  the  final  distribution  of  p',  o',  T'  and  u1  in 
the  JPL/HWT  boundary- layer  is  shown  on  Figure  3.  Generally  the  fluctuations  peak 
at  the  sublayer  edge  as  well  as  in  the  main  boundary- lay er  itself.  Comments  on 
extrapolating  or  interpolating  these  results  to  different  flow  conditions  are 
given  in  Section  IV. 


Some  other  recent  findings  are  also  important.  One,  Bhown  on  Figure  4,  is  that 
the  turbulence  Reynolds  number  Re^  is  extremely  high.  Since  a  high  Re^  denotes 
also  a  highly  equilibrated  layer,  and  since  the  suspicions  about  the  pressure 
gradient  3p/3y  have  now  been  allayed  (Section  2.2.1),  the  extreme  "maturity"  of 
this  boundary-layer  has  baen  confirmed. 


Figure  5  plots  the  final  picture  of  the  distribution  of  integral  lengths  A. 
Repetitive  tests  with  different  sensors,  interpreted  in  different  ways,  are  in 
mutual  agreement  that  Au  and  are  concentrated  at  A/6  ■  1.8/10  ■  0,18  with  a 

"scatter"  band  of  about  30*/..  However,  the  pressure  scale  Ap  is  much  smaller. 
Presently  the  magnitude  of  Ap  shown  on  Figure.  5  is  somewhat  unsettled  for  the 
following  reason.  Like  Au  and  Aft  Ap  is  computed  from  the  data  via  Taylor's  hy¬ 
pothesis,  i.e., 


A 


P 


ii  fp--2-Cf  -  0)s 
4  1  .2  ; 


2 

where  u  is  the  local  velocity,  p'  (f  ■  0)  the  low-frequency  spectral  limit  and 
p'  the  wideband  pressure  fluctuations.  Within  the  test  accuracy,  pl2(f  «  0)/p'2 
(which  is  an  integral  time  scale)  is  definitely  known,  but  the  factur  u/4  used  to 
convert  to  a  space  scale  is,  perhaps,  used  too  prematurely.  The  point  is  that 
the  propagation  speed  of  p'  is  radically  different  from  that  of  u'  and  T1,  which 
are  known  to  convect  with  the  local  flow  speed  (Reference  4),  Therefore  the 
question  of  Ap  cannot  be  finally  settled  unless  a  suitable  propagation  velocity 
is  found  to  replace  u. 


2.3  STUDY  OF  TURBULENCE  AT  REDUCED  REYNOLDS  NUMBER 

In  addition  to  the  TBL-A  results  presented  above  (and  in  Reference  1)  which  were 
done  at  PQ  -  3200  cm  Hg  (Re  «  49,500  per  cm),  some  data  are  available  on  tape 
from  tests  at  PQ  -  1700  cm  Hg  (Re  ■  26,300  per  cm).  In  the  present  contract  period 
these  data,  too,  were  analysed. 

The  experimental  procedures,  instrumentation  and  method  of  data  reduction  were 
identical  to  those  described  in  Reference  1.  The  measurements  were  made  using 
hot-wire  6-3,  the  Bame  anemometer  with  which  approximately  one-half  the  turbu¬ 
lence  data  at  P0  -  3200  cm  Hg  was  acquired.  Calculation  of  fluctuation  intensities 
included  the  effect  of  non-zero  pressure  fluctuations  using  the  assumption  Rai|  ■ 

R.rn  ■  0,  RoT  ■  -1.  The  wideband  fluctuation  intensities  of  velocity,  temperature 
and  pressure,  which  were  obtained  from  integrating  over  frequency  the  results  of 
spectral  modal  analysis,  are  listed  below  as  a  function  of  distance  from  the  wall, 
y.  At  each  y  position  both  the  spectral  results  and  the  wideband  mode  diagram 
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FIGURE  3.  VARIATION  OF  PRESSURE,  DENSITY,  TEMPERATURE  AND  VELOCITY 
FLUCTUATION  INTENSITIES  IN  THE  BOUNDARY  LAYER. 


FIGURE  4.  TURBULENCE  REYNOLDS  NUMBER  VARIATION  ACROSS  THE  HYPERSONIC 
BOUNDARY  LAYER. 


LENGTH  distribution. 


were  carefully  examined  to  assess  the  data  and  three  points  were  found  to  be  un¬ 
acceptable.  At  y  *  0.968  Inches  an  excessive  scatter  in  the  data  was  observed 
which  precluded  a  rational  curve  fit  to  find  the  fluctuation  Intensities.  While 
the  data  was  better  behaved  at  y  *  4.469  inches,  errors  in  several  low  current 
data  points  resulted  in  negative  mean  square  values  for  the  turbulent  fluctua¬ 
tions.  Although  the  results  of  this  point  could  be  improved  by  eliminating  the 
offending  low  current  data  points,  this  procedure  was  not  carried  out.  Finally, 
at  y  ■  2.967  inches,  an  error  in  the  lowest  current  data  point  led  to  an  over¬ 
prediction  in  the  pressure  fluctuations  by  50%.  Since  this  y  position  is  located 
near  the  point  of  peak  fluctuations,  a  wideband  analyais  was  carried  out  with  the 
erroneous  point  omitted,  yielding  the  corrected  results  listed  below. 


P  ■  1700  cm  WIDEBAND  FLUCTUATIONS 
o 


1 

H 

T1  (%) 

2im 

0.96  b 

0.92 

2.52 

5.49 

1.467 

0.98 

3.76 

7.80 

1.968 

0.80 

3.92 

8.35 

2.464 

0.72 

4.35 

9.79 

2.967 

0.61 

4.70 

12.00 

3.467 

0.45 

4.88 

11.53 

3.969 

0.25 

3.26 

7.64 

4.469 

0,07 

0.96 

2.88 

4.968 

0.031 

1.23 

2.99 

The  seven  acceptable  results  for  velocity,  temperature  and  pressure  fluctuations 
are  plotted  in  Figure  6  for  comparison  with  the  3200  cm  Hg  results,  which  are 
represented  by  curves  re-drawn  from  Figure  3.  It  is  seen  that  near  the  outer 
edge  of  the  boundary- layer  (y  a>  3.00")  the  resultB  at  the  two  pressures  are  in 
good  agreement.  Closer  to  the  wall,  however,  the  1700  cm  Hg  results  yield  fluc¬ 
tuation  intensities  that  are  less  than  those  at  3200  cm  Hg,  particularly  the 
temperature  fluctuations  which  differ  by  as  much  as  a  factor  of  two.  In  making 
such  a  comparison  it  should  be  recalled  that  of  the  two  hot-wires  used  at  3200 
cm  Hg,  wire  6-3  consistently  gave  results  that  were  slightly  smaller  and  this 
may  be  reflected  in  the  results  at  1700  cm  Hg. 

Figures  7  and  8  show,  respectively,  the  variation  of  the  quantity  T'/(Y-l)M2u ' 
and  the  integral  scale  lengths,  A,  as  a  function  of  y.  In  foth  figures,  the 
agreement  between  the  results  obtained  at  the  two  tunnel  pressures  is  quite  good. 

In  view  of  this  it  is  not  clear  whether  the  difference  in  fluctuation  intensities 
obtained  within  the  interior  of  the  boundary-layer  is  due  to  the  change  in  Reynolds 
number  corresponding  to  the  two  tunnel  pressures,  or  is  a  consequence  of  the  pecul- 
arities  of  the  hot-wire  with  which  the  data  was  acquired.  Since  the  former  alter¬ 
native  is  in  variance  with  other  evidence  (see  Section  IV)  it  raises  a  need  to  con¬ 
tinue  and  intensify  the  search  for  the  leynolds  number  effect  on  the  boundary- 
layer  turbulence. 


FIGURE  6.  FLUCTUATION  INTENSITY  DISTRIBUTION  IN  THE  JFL/3WT  BOl 
LAYER  AT  P_  -=  1700  cm  Hg  and  3200  cm  Eg- 


FIGURE  7.  COMPARISON  OF  THE  STRONG  REYNOLDS  ANALOGY  IN  THE  BOl 
LAYER  FOR  1700  cm  Hg  AND  3200  cm  Hg. 


FIGURE  8.  COMPARISON  OF  SCALE  LENGTHS  FOR  THE  1700  cm  AND  3200 


SECTION  III 


HYPERSONIC  CONE  BOUNDARY- LAYER  EXPERIMENT 


3.1  OBJECTIVE 

The  purpose  of  this  experiment  is  to  determine  the  mean  and  fluctuating  charact¬ 
eristics  of  a  hypersonic  turbulent  boundary-layer  on  a  slender  porous  cone,  The 
tost  was  conducted  at  the  Arnold  Engineering  and  Development  Center  (AEDC) ,  Von 
Karman  Facility,  using  Tunnel-B  at  Mach  8,  This  experiment  represents  a  joint 
effort  between  the  General  Electric  Re-Entry  and  Environmental  Systems  Division 
and  the  Aeronutronic  Division  of  Philco-Ford,  originally  with  G.E.  responsible 
for  the  mean-flow  measurements  and  Philco-Ford  for  the  fluctuation-flow  measure¬ 
ments  . 

3.2  HOT-WIRE  SURVIVAL  TESTS 

Since  the  AEDC  Tunnel  B  represented  a  now  environment  for  the  Philco-Ford  designed 
hot-wire  probea,  and  because  the  flow  conditions  are  more  severe  than  those  en¬ 
countered  in  both  the  Philco-Ford  SWT  and  the  JPL/HWT,  a  series  of  testB  to  examine 
the  survivability  and  endurance  of  the  probes  was  carried  out  over  the  eight-month 
period  prior  to  the  first  Philco-Ford  entry.  To  accelerate  the  evaluation  process, 
an  instrumentation  rake  with  provisions  for  mounting  8  probes  was  used,  with  each 
probe  equipped  with  an  individual  current  monitoring  circuit,  The  output  from 
each  circuit  was  fed  to  a  strip  chart  recorder  in  order  to  provide  a  complete  time 
history  of  the  probe  during  its  exposure  to  the  flow.  Details  of  the  rake  and  itu 
operation  are  described  in  Reference  5.  A  photograph  of  the  rake  located  in  the 
free-etream  flow  of  Tunnel  B  is  shown  in  Figure  9. 

An  initial  series  of  tests  was  carried  out  on  a  parasitic  basis  with  the  coopera¬ 
tion  and  assistance  of  ARO  personnel  during  the  period  from  January  -  March  1972. 

On  the  basis  of  these  initial  tests,  which  are  described  in  Reference  5,  It  was 
concluded  that: 

(1)  Survival  time  of  the  hot-wires,  particularly  the  0.000010  inch  dia¬ 
meter  probes,  was  too  short  to  be  of  practical  value  for  acquisition 
of  turbulence  data, 

(2)  Survivability  improves  with  Increasing  wire  diameter  and  decreasing 
aspect  ratio, 

(3)  Shielding  of  the  hot-wire  is  essential  during  tunnel  start  and  dur¬ 
ing  changes  in  stagnation  pressure. 

(4)  The  stagnation  temperature  muat  be  maintained  above  the  condensation 
threshold  when  the  wire  is  expoued  to  the  flow. 

Based  on  these  and  subsequent  test  results  it  was  decided  to  qualify  all  hot¬ 
wires  in  the  Philco-Ford  SWT  prior  to  use  in  Tunnel  B  in  order  to  eliminate 
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FIGURE  9.  INSTRUMENTATION  Rj 
LOCATED  IB  FREE-S' 


structurally  weak  probes.  In  addition,  during  the  qualifying  cast,  the  probe 
output  was  monitored  and  those  which  exhibited  erratic  or  abrupt  variations  in 
output  signal  were  discarded. 

A  second  series  of  survivability  tests  was  carried  out  in  May  1972  during  the 
Phase  I  entry  of  the  G. E/Phi lco-Ford  experiment.  The  Phase  I  test,  represent¬ 
ing  the  first  installation  of  the  cone  model  in  Tunnel  B,  was  conducted  by  G.E. 
to  check-out  the  test  hardware  and  obtAin  preliminary  mean-flow  measurements. 

On  the  first  hot-wire  test,  eight  probes  consisting  of  four  10  pinch,  two  20 
pinch  and  two  50  pinch  diameter  wires  were  installed  in  the  tunnel,  and  tested 
in  the  free-strcam  flow  at  supply  pressures,  PQ,  of  200,  400,  600,  and  850  paia. 
The  supply  pressure  was  initially  established  at  200  psla  and  after  a  suitable 
interval  increased  approximately  200  psia  until  a  P  of  850  psia  was  reached. 

At  each  P0  the  supply  temperature  was  maintained  sufficiently  high  (760  -  890°F) 
to  avoid  condensation  in  the  flow  and  during  tunnel  start  and  change  in  supply 
conditions  the  probes  were  retracted  into  the  recess  in  the  ceiling  of  the  tunnel 
for  protection.  At  a  supply  pressure  of  200  and  400  psia  all  the  wires  survived 
exposure  to  the  flow  for  a  duration  of  approximately  45  minutes  at  each  condi¬ 
tion.  With  P0  ■  600  psia,  two  10  iiinch  and  both  20  (iinch  diameter  wires  failed 
after  10-17  minutes,  while  the  remaining  wires  survived  for  64  minutes.  The 
supply  pressure  was  then  increased  to  850  psia,  where  the  remaining  10  iiinch 
wires  failed  after  4  to  19  minutes,  while  the  50  (iinch  wires  survived  for  ap¬ 
proximately  one  hour  before  failure  occurred. 

On  a  second  hot-wire  test,  four  10  ginch  and  four  20  iiinch  diameter  wires  were 
installed  in  the  tunnel  and  the  test  was  conducted  at  a  P0  of  850  psia.  While 
the  10  iiinch  wires  failed  after  durations  hanging  from  3  to  9  minutes,  greater 
success  was  obtained  with  the  20  (iinch  wires,  with  one  probe  surviving  for  as 
long  as  100  minutes.  During  this  entry,  heat  transfer  tests  of  the  G.E.  cone 
model  were  carried  out  and  it  was  discovered  that  the  model  could  be  inserted 
into  and  retracted  from  the  tunnel  without  causing  either  mechanical  or  elec¬ 
trical  damage  to  the  hot-wire  probes.  This  finding  was  of  practical  importance 
for  the  conduct  of  the  final  tests  since  it  minimized  the  precautions  needed  to 
protect  the  hot-wire  probe  during  insertion  of  the  cone  model. 

A  final  survivability  test,  conducted  on  a  parasitic  basis  one  week  after  the 
Phase  I  entry  used  four  10  (iinch  and  four  50  (iinch  diameter  wires.  Two  of  the 
small  wires  were  lost  prior  to  the  tunnel  start,  apparently  due  to  burn-out 
caused  by  pickup  of  a  transient  voltage  pulse,  and  one  50  (iinch  wire  was  lost 
when  inserted  into  the  flow.  The  remaining  wires  were  successfully  inserted 
into  the  flow  on  two  occasions  and  survived  a  total  of  70  minutes  at  tunnel 
supply  conditions  of  800  psia  and  875°F.  Although  these  wires  were  intact  at 
the  end  of  the  test  and  survived  retraction  into  the  top  recess  of  the  tunnel, 
the  two  10  (iinch  wires  failed  during  or  immediately  following  tunnel  shut-down. 


*  Initial  difficulties  are  best  understood  when  one  considers  that  the  least 
delicate  wires  listed  here,  of  50  (iinch  diameter,  generally  represented  the 
forefront  of  the  state-of-the-art  at  the  beginning  of  this  work. 
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On  the  basis  of  chase  tests, summarized  in  Table  I,  the  10  p,inch  diameter  wires, 
as  well  as  the  20  and  50  p,inch  wires,  demonstrated  excellent  survivability  and 
endurance  at  a  tunnel  supply  pressure  of  400  psia,  the  lowest  PQ  of  the  final 
test  matrix.  On  the  other  hand,  at  850  psia,  the  highest  P  to  be  used  in  the 
final  tests,  only  the  50  p.inch  wires  showed  acceptable  survivability  and  endur¬ 
ance  while  the  smaller  wires  Indicated  improved  but  still  questionable  perfor¬ 
mance.  It  was  concluded,  therefore,  that  the  final  turbulence  measurements 
could  be  carried  out  using  wires  as  small  as  10  p,  at  the  400  psia  condition  and 
with  50  and  possibly  20  pinch  wires  at  850  psia.  In  all  cases  it  was  considered 
desirable  to  restrict  the  aspect  ratio  of  the  wire  to  less  than  200. 

A  final  series  of  parasite  tests  of  10,  20  and  50  pinch  diameter  hot-wires  was 
conducted  in  AEDC  Tunnel  B  during  the  last  week  in  July  1972.  The  tunnel  was 
operated  at  Mach  8  with  stagnation  conditions  of  800  psia  and  875°F  and  with  the 
probes  located  in  the  free-stream.  The  purpose  of  these  tests  was  to  determine 
experimentally  the  overheat  current  at  which  hot-wire  failure  occurred  as  a  func¬ 
tion  of  hot-wire  diameter.  The  burn-out  currents  were  converted  to  wire  temper¬ 
atures  and  indicated  that  burn-out  occurred  when  the  hot-wire  was  heated  to  1000- 
],100°K.  This  result  is  slightly  less  than  that  observed  during  the  JP.L  tests  al¬ 
though  the  dynamic  pressure,  and  consequently  the  structural  load  on  the  hot-wire, 
is  larger  in  the  AEDC  facility.  Consequently  the  observed  burn-out  currents  were 
consistent  with  the  JPL  findings. 

During  this  last  test  the  hot-wire  rake  was  modified  to  accommodate  four  total 
temperature  probes  of  different  construction  which  had  been  designed  by  G.E. 
Previous  designs  used  by  G.E.  had  been  found  to  have  recover  factors  signifi¬ 
cantly  less  than  unity  which  increased  the  number  of  iterations  required  to  find 
the  true  total  temperature.  The  new  T0  probes  were  tested,  therefore,  to  deter¬ 
mine  the  effect  of  probe  geometry  on  its  recovery  factor  characteristic.  The 
tests  were  successful  in  clearly  identifying  a  probe  with  superior  character¬ 
istics,  and  on  this  basis,  a  T0  probe  to  be  used  by  both  Philco-Ford  and  G.E. 
during  the  final  tests  was  selected, 


3.3  FINAL  EXPERIMENT  DESIGN 
3.3.1  TEST  FACILITY 

The  present  investigation  was  conducted  at  the  Arnold  Engineering  Development 
Center  (AEDC),  Von  Karman  Facility,  Tunnel  B  at  Mach  8.  This  facility  is  a  con¬ 
tinuous  flow,  hypersonic  wind  tunnel  with  a  50-inch  diameter  test  section.  Inter¬ 
changeable  axisyinmetric  contoured  nozzles  provided  flow  at  either  Mach  6  or  8.  The 
tunnel  is  a  closed-circuit  type  which  operates  with  air  supplied  by  a  central  com¬ 
pressor  system  over  stagnation  pressures  ranging  from  50  to  850  psia.  The  corres¬ 
ponding  flow  conditions  provide  free-stream  Reynolds  numbers  from  0.30  x  10(l  to 
3.8  x  10^  per  foot  at  Mach  8  with  a  total  temperature  of  1350°F,  The  entire  tunnel 
is  cooled  by  integral,  external  '  ter  Jackets.  Complete  specification  he 

facility  are  described  in  Reference  6,  while  the  tunnel  is  shown  schei  ly  in 

Figure  10. 


TABLE  I 


RESULTS  OF  HOT-WIRE  ENDURANCE  TESTS  -  MAY  18-23 
Log  of  Time  (in  minutes)  Exposed  to  Flow 


WIRE 

NUMBER 

WIRE  DIAMETER 
(finches) 

ASPECT 

RATIO 

SUPPLY  PRESSURE, 
200  400  600 

P0,  PSIA 
850(a) 

TOTAL 

TIME 

23 

50 

220 

45 

48 

64 

60* 

217 

22 

50 

210 

45 

48 

64 

63* 

210 

15 

20 

190 

45 

48 

17* 

110 

19 

20 

290 

45 

48 

10* 

103 

9 

10 

180 

45 

48 

64 

19* 

176 

4 

10 

260 

45 

48 

17* 

110 

6 

10 

280 

45 

48 

64 

4* 

161 

7 

10 

420 

45 

48 

14* 

107 

3 

10 

265 

3* 

3 

21 

10 

245 

9* 

9 

10 

10 

220 

6* 

6 

30 

10 

260 

7* 

7 

17 

20 

215 

4* 

4 

16 

20 

190 

20* 

20 

28 

20 

260 

58* 

58 

29 

20 

205 

100* 

100 

11 

10 

285 

+ 

14 

10 

255 

+ 

■  13 

10 

175 

70 

70** 

31 

10 

350 

70 

70** 

24 

50 

180 

70 

70** 

25 

50 

240 

1 

26 

50 

165 

70 

70** 

27 

50 

155 

70 

70** 

*  Denotes  time  of  failure  after 

exposur 

e  at 

indicated  P  . 

o 

**  Wire 

intact  at  end 

of  test. 

\ 

Lost 

during 

insertion 

into  flow. 

+  lost 

prior  to  tunnel  start. 

(a) 

P  - 

0 

800  for  last  8  wires. 

FIGURE  10.  SCHEMATIC  OF  AEDC  TUNNEL  B 


In  order  to  provide  quick  access  to  the  model  without  interrupting  the  tunnel 
operation,  a  vacuum  sealed  tank  is  located  immediately  below  the  test  section 
Figure  10,  into  which  the  model  and  its  support  can  be  retracted.  When  the  model 
is  retracted  the  test  section  is  sealed  from  the  tank  which  can  then  be  vented  to 
atmosphere.  After  the  tank  entrance  door  is  closed,  the  tank  is  vented  to  the 
test  section  pressure,  the  door  between  the  tank  and  section  is  opened,  and  the 
model  is  injected  into  the  airstream. 

Because  of  a  similar  need  to  provide  quick  access  to  the  hot-wire  probe,  it  was 
necessary,  prior  to  the  fluctuation  measurements,  to  construct  a  vacuum  sealed 
tank  into  which  the  probe  traversing  mechanism  could  be  retracted.  This  tank 
was  located  above  the  test  section  over  an  existing  port  used  for  observation 
and  photography  of  the  model  (see  Figure  10).  The  tank,  shown  in  its  operating 
position  in  Figure  11a,  is  raised  approximately  6  feet  above  the  top  of  the  test 
section  by  means  of  an  overhead  hoist  to  permit  access  to  the  probe  and  its 
traversing  and  retraction  mechanism,  Figure  lib. 

3.3.2  TEST  HARDWARE 

The  model  used  in  this  investigation,  which  was  supplied  by  G.E.,  is  a  3°  half¬ 
angle  cone,  60  inches  long  (10.46  inch  base  diameter).  It  consists  of  an  imper¬ 
vious  sharp  nose  which  is  10  Inches  long,  and  a  porous  frustum  which  is  50  inches 
long.  The  porous  section  of  the  model  is  subdivided  into  four  independent  com¬ 
partments,  each  with  separate  gas  supply  tubes  which  are  1/4"  O.D.  stainless  steel 
tubing.  The  forward  two  chambers  each  have  one  supply  tube,  the  rear  two  chambers 
each  have  two  supply  tubes.  A  photograph  qf  the  model  assembly  located  in  the 
tank  below>  the  test  section  is  shown  in  Figure  12. 

The  porous  frustums  of  the  model,  which  were  fabricated  at  Mott  Metallurgical 
Corp.,  Farmington,  Conn.,  are  made  of  316  sintered  stainless  steel  which  is  ap¬ 
proximately  60*4  dense.  The  wall  is  approximately  1/4"  thick.  These  porous  ele¬ 
ments  are  supported  by  a  stainless  steel  substructure  which  contains  all  of  the 
necessary  internal  plumbing  and  instrumentation  routing  ports.  The  cone  is  in¬ 
strumented  with  26  heat  transfer  gages  and  34  pressure  orifices,  the  locations  of 
which  are  indicated  in  Figure  13.  The  surface  heat  transfer  measurements  were 
made  with  a  copper-constantan  Gardon  type  heat,  gage  which  provided  an  indication 
of  both  the  heat  transfer  rate  and  the  local  surface  temperature, 

Calibration  of  both  the  gross  porosity  of  each  chamber  of  the  model  and  of  the 
local  mass  flux  variation  have  been  made  at  G.E.  prior  to  the  test  program. 

These  measurements  indicate  that  the  local  variations  do  not  exceed  the  ±  107, 
design  specification,  In  addition,  the  grosB  porosity  of  each  chamber:  was  veri¬ 
fied  at  AEDC  before  and  after  each  entry,  thereby  indicating  that  the  porosity 
of  the  chambers  remained  constant  during  the  test.  The  gas  metering  systems 
provided  by  AEDC  for  each  chamber  were  accurate  to  within  several  percent  of 
the  desired  test  condition. 

A  detailed  description  of  the  model  and  the  surface  instrumentation  can  be  found 
in  Reference  7. 
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FIGURE  12.  PHOTOGRAPH  OF  GE  CONE  MODEL  LOCATED  IN  ACCESS  TANK 


FIGURE  13.  LOCATION  OF  CONE  SURFACE  PRESSURE  AND  TEMPERATURE  GAUGES. 


3,3.3  INSTRUMENTATION 


Although  G.E.  was  responsible  for  characterizing  the  mean  boundary- layer  flow, 
it  was  decided  to  make  and  tape-record  mean-flow  measurements  during  the  Philco- 
Ford  entry  in  order  to  provide: 

(1)  A  check  on  the  more  detailed  G.E.  measurements,  thereby  verifying 
the  constancy  of  the  boundary- layer  flow,  and 

(2)  A  complete  permanent  set  of  boundary-layer  data  at  all  conditions 
of  the  Philco-Ford  test  matrix. 

The  latter  was  particularly  crucial  since,  in  spite  of  initial  intentions  to 
carry  out  the  mean  and  fluctuating  flow  measurements  at  Identical  teBt  conditions, 
the  final  G.E,  and  Philco-Ford  tests  were  in  fact  performed  under  different  con¬ 
ditions  . 


a.  Instrumentation  Rake.  The  Philco-Ford  instrumentation  rake,  shown  in 
Figure  14,  consisted  of  four  1/8  inch  O.D,  x  3  Inch  long  stainlcsB  steel  tubes 
aligned  in  a  row  and  silver  soldered  to  each  other  and  to  a  single  wedge  shaped 
sting  which  mounted  to  the  AEDC  probe  traversing  mechanism.  The  probes  were  in¬ 
stalled  on  the  rake  by  insertion  through  the  1/8  inch  diameter  tubes,  with  a 
pitot  probe  and  total  temperature  probe  occupying  the  outboard  positions  and  a 
hot-wire  probe  and  surface  sense  probe  located  in  the  inboard  positions.  The 
probes  were  aligned  with  their  sensing  tips  essentially  the  same  distance  from 
the  cone  tip,  but  with  the  pitot,  total  temperature  and  hot-wire  probes  dis¬ 
placed  from  0.005  to  0.040  inches  vertically  above  the  surface  sensor  in  order 
to  avoid  contact  with  the  cone  surface.  The  lateral  and  vertical  displacement 
of  each  probe  (relative  to  the  surface  sense  probe)  and  the  effect  of  surface 
curvature  were  taken  into  account  to  convert  the  indicated  probe  position  into 
its  true  height  above  the  surface  along  a  line  passing  through  the  center  of  the 
cone.  This  procedure  is  described  in  Appendix  A. 

b.  Probe  Design.  The  hot-wire  anemometer  consisted  of  a  Pt-107.  Rh  wire 
mounted  across  the  tips  of  two  sharp-pointed  prongB  and  bonded  to  the  supports 
with  a  commercial  gold  paint.  For  the  present  experiments  wires  of  10,  20  and 
50  |j,inch  diameter  were  used  and  during  the  week  prior  to  the  test  approximately 
one  dozen  probes  of  each  size  were  fabricated.  The  wires  were  mounted  with  suf¬ 
ficient  slack  to  eliminate  spurious  signals  due  to  "strain  gage"  effects  and  the 
wire  aspect  ratio  was  restricted  to  the  range  from  150  to  250. 

The  probe  body  consisted  of  a  3/32  inch  diameter  x  1-1/4  inch  long  alumina  cylin¬ 
der  cemented  at  one  end  to  a  1/32  inch  diameter  x  1/4  inch  long  alumina  cylinder 
uaed  to  support  the  hot-wire  prongs.  Both  cylinders  were  provided  with  two  holes 
along  their  length  through  which  the  hot-wire  leads  were  inserted.  The  leads 
were  30  gage  copper  wire  encased  in  high  temperature  fiberglass  insulation.  Three- 
mil  nickel  wire,  used  for  the  hot-wire  prongs,  was  silver-soldered  to  one  end  of 
the  vopj'-tr  leads.  The  prongs,  leads,  and  alumina  tubes  were  held  together  with 
i  *»>;  vktss  Company  Pyroceram  cement.  Details  of  the  probe  construction  are 
•  in  Reference  8.  A  photograph  of  the  assembled  probe  is  shown  in  Figure 
!  ,  c  an  enlargement  of  the  probe  tip,  showing  a  20  uinch  diameter  wire 

rth'J:U  rl.  U'«i  the  prongs,  is  presented  in  Figure  15b. 
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FIGURE  14.  PHI LCO- FOR D  INSTRUMENTATION  RAKE  WITH  PROBES  INSTALLED. 

FROM  FRONT  TO  BACK:  PITOT  PROBE,  HOT-WIRE,  SURFACE  SENSOR 
TOTAL  TEMPERATURE  PROBE. 


After  the  probe  was  assembled  and  qualified  in  the  Philco-Ford  SWT,  the  hot-wire 
was  oven  calibrated  to  determine  its  temperatur j-resistance  characteristics. 
Following  the  oven  calibration,  the  wire  resistance  was  checked  daily  prior  to 
its  use  in  Tunnel  B  to  ensure  that  the  wire  characteristics  remained  unchanged. 

The  total  temperature  probe  consisted  of  a  chromel-nlumel  thermocouple,  sheathed 
within  q  0.020  inch  O.D.  stainless  steel  tube,  which  was  inserted  within  a  0.030 
inch  0.1).  tube  for  rigidity.  The  0.030  inch  tube  in  turn  was  silver  soldered 
within  a  C.U93  inch  O.D.  tube  for  insertion  into  the  rake.  The  sensing  end  of  the 
probe  was  shielded  by  a  0.030  inch  O.D.  tube  as  shown  in  Figure  16. 

The  pitot  probe  was  constructed  in  a  similar  fashion  with  a  0.030  inch  O.D.  tube 
transitioning  to  a  0.093  inch  O.D.  tube  for  installation  in  the  instrumentation 
rake,  The  sensing  end  of  this  probe  was  flattened  to  an  outside  height  of  0,020 
inches  as  indicated  in  Figure  16. 

Finally,  the  surface  sense  probe,  also  shown  in  Figure  16,  consisted  of  a  0.015 

inch  diameter  hardened  steel  wire  held  in  place  with  Sauereiaan  cement, to  provide 

electrical  isolation, within  a  0.093  inch  O.D,  tube.  The  sense  probe  was  connected 
in  series  with  a  low  voltage  battery  and  a  lamp,  which  was  turned  on  when  the 
probe  contacted  the  cone  to  complete  the  circuit. 

c.  Probe  Traverse  Mechanism.  The  probe  traversing  mechanism  constructed 
for  the  fluctuating  flow  measurements  wsb  designed  to  permit  retraction  out  of  the 
flow  and  allow  access  to  the  probeo  without  interrupting  the  tunnel  operation.  In 
addition,  the  traverse  mechanism  provided  for  remote  drive  in  the  direction  parallel 
to  the  tunnel  axis  (x)  and  in  the  direction  normal  to  the  cone  surface  (y).  The 
axial  location  of  the  probes  was  defined  by  locating  the  sense  probe  adjacent  to 
the  cone  at  the  Junction  between  two  chambers  to  establish  an  x  reference  immed¬ 
iately  after  inserting  the  probes  into  the  flow.  The  probe  could  then  be  located 

at  a  desired  x  station  by  driving  the  traverse  mechanism  in  the  required  distance* 

from  the  known  reference  position.  The  axial  station  could  be  determined  within  an 
accuracy  of  ±  0.05  inches  and  a  repeatability  of  4  0.10  Inches.  The  vertical  height 
above  the  model  was  determined  by  driving  the  probe  toward  the  cone  surface  and 
establishing  a  y  -  0  position  with  the  surface  sense  probe.  All  boundary- layer 
surveys  were  made  driving  away  from  the  surface  after  the  y  ■  0  location  was  iden¬ 
tified.  The  y  station  could  be  determined  within  4  0.001  inches  and  was  repeat- 
able  within  4  0.005  inches. 

In  order  to  protect  the  hot-wire  probe  during  insertion  into  the*  tunnel,  change  of 
tunnel  conditions,  etc.,  the  traverse  mechanism  was  provided  with  a  pneumatically 
actuated  shield  which  covered  the  entire  front  portion  of  the  instrumentation  rake, 
Figure  17  shows  the  probe  traverse  mechanism  located  atop  the  tunnel  test  section 
with  the  shield  in  both  its  lowered  and  raised  positions,  Figure  18  shows  the 
probes  located  within  the  cone  boundary- layer  and  indicates  the  clearance  between 
the  cone  and  the  shield  when  the  latter  is  in  its  retracted  position. 

3.3.4  TEST  PROCEDURE 

a.  Test  Matrix.  Several  discrete  entries  in  Tunnel  11  were  required  for  com¬ 
pletion  of  both  the  mean  and  fluctuating  flow  measurements.  The  measurements  were 
carried  out  in  accordance  with  the  following  schedule: 
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FIGURE  17.  PHOTOGRAPH  OF  PROBE  TRAVERSING  MECHANISM  LOCATED  OR  TOP 
OF  TEST  SECTION  SHOWING  HOT-WIRE  SHIELD  IN  RAISED  AND 
LOWERED  POSITIONS. 


ENTRY 


PURPOSE 


Phase  I 

Phase  II 
Phase  III 

Phase  IV 


Heat  transfer  measurements, 
hardware  check-out  and  hot-wire 
survivability . 

Mean  flow  profiles, 

Fluctuating  flow  and  mean  flow 
profiles . 

Fluctuating  flow  and  mean  flow 
profiles. 


Originally  the  turbulence  measurements  were  to  be  carried  out  for  the  conditions 
listed  in  Table  II,  which  corresponded  to  the  baseline  conditions  of  the  G.E. 
mean  flow  measurements.  However,  because  of  difficulties  encountered  with  the 
hot-wire  at  BSO  psia,  this  condition  was  relaxed  and  final  tests  performed  at 
600  psia.  The  completed  test  matrix  is  summarized  in  Table  III.  Since  the  G.E. 
program  was  carried  out  only  at  850  psia  (omitting  the  originally  planned  350 
psia  tests),  the  mean  flow  data  acquired  by  Fhllto-Ford  during  the  hot-wire  tests 
represents  a  necessary  ingredient  for  the  successful  completion  of  thiB  task. 


b.  Measurement  Techniques.  Although  the  mean  flow  measurements,  probe 
positions,  and  tunnel  conditions  were  recorded  and  processed  by  the  AEDC  Data 
Acquisition  System,  the  probe  data  was  also  tape-recorded  by  Philco-Ford  in  order 
to  provide,  with  the  hot-wire  measurements,  a  complete  magnetic  tape  data  file. 
Block  diagrams  of  the  Philco-Ford  data  recording  schemes  are  shown  in  Figures  19 
and  20.  These  schemes,  called  configurations  A  and  B,  differ  only  in  the  manner 
in  which  the  hot-wire  signal  was  processed.  In  both  cases  the  pitot  pressure, 
total  temperature  and  surface  sense  probe  positions  were  monitored  by  DVM's  and 
converted  to  frequency  signals  for  tape  recording  by  the  Honeywell  7620  recorder. 
The  hot-wire  current  and  mean  voltage  were  also  monitored  by  DVM's  and  the  latter 
converted  to  a  frequency  signal  for  tape  recording. 


The  fluctuating  voltage  output  from  the  hot-wire  was  monitored  by  a  Hewlett-Packard 
wave  analyzer  an^)  Tektronix  Oscilloscope  in  order  to  observe  the  frequency  content 
and  quality  of  tone  signal.  With  configuration  A  the  signal  was  then  directly  re¬ 
corded,  while  vjith  configuration  B  the  signal  was  fed  to  An  RMS  muter  and  converted 
into  a .wideband,  signal  prior  to  tape  recording.*  The  purpose  of  the  two  recording 
schemes  is  described  below. 


*  The  hot-wire  amplifier  was  operated  typically  at  a  high  frequency  cut-off  of 
500  kHz,  Since  the  FM  channel  cuts  off  at  400  kHz  when  the  recorder  is  oper¬ 
ated  at  120  ipa,  the  signal  was  also  recorded  on  a  direct  channel  which  extends 
to  1.5  MHz.  Thus,  although,  the  direct  recording  is  of  poorer  quality  than  the 
FM  recording,  this  procedure  allows  the  complete  amplifier  output  to  be  recorded 


ORIGINAL  TEST  MATRIX  FOR  TURBULENCE  MEASI 


final  test  matrix  for  turb 


FIGURE  19.  BLOCK  DIAGRAM  OF  TURBULENCE  HEAS 
IN  CONFIGURATION  A. 


CONFIGURATION  B. 


Prior  tc,  the  Phase  III  tests,  it  was  decided  to  cover  that  part  of  the  test 
matrix  which  ensured  maximum  survivability  of  the  hot-wires.  For  this  reason, 
the  "safest"  combination  of  wire  diameters  (O.OOOQ5inches)  and  tunnel  pressures 
(P0  ■  350  psia)  were  used.  For  this  same  reason,  it  was  also  decided  to  init¬ 
ially  operate  the  wires  using  configuration  B.  This  configuration  involves 
traversing  the  hot-wire  across  the  boundary-layer  continuously  at  constant  cur¬ 
rent  and  recording  at  7^  ips,  so  that  the  wideband  fluctuation  intensity  distri¬ 
bution  with  vertical  distance  could  he  extracted  from  the  data.  With  config¬ 
uration  A,  the  preferred  recording  scheme,  the  hot-wire  is  located  at  several 
pre-selected  distances  above  the  surface  and  its  output  recorded  at  120  ips  as 
the  heating  current  is  changed  in  discrete  steps.  The  latter  has  the  advantage 
over  configuration  B  in  that  it  can  yield  spectral  information  and  thus  more 
accurate  fluctuation  data;  but,  on  the  other  hand,  it  Involves  greater  risk 
than  configuration  B  because  it  takes  data  at  a  slower  rate. 

This  approach  was  initiated  during  the  Phase  III  tests.  It  was  further  planned 
that  when  data  acquisition  using  configuration  B  was  completed  for  a  given  test 
condition,  the  test  would  be  repeated  using  configuration  A  in  order  to  provide 
the  more  accurate  fluctuation  data.  However,  for  reasons  discussed  later,  little 
quantity  information  was  collected  during  the  Phase  III  entry.  In  addition,  on 
the  basis  of  the  experience  acquired  with  the  hot-wires  during  this  entry,  it  was 
believed  that  the  wires  would  survive  for  significant  periods  of  time.  There¬ 
fore,  configuration  B  was  deleted  for  the  Phase  IV  tests,  and  the  test  matrix 
outlined  in  Table  III  was  completed  using  configuration  A.  Figure  21  shows  the 
layout  of  the  Philco-Ford  electronics  located  just  downstream  of  the  tunnel  test 
section. 

3.4  TEST  RESULTS 
3.4.1  PHASE  III  TESTS 

A  summary  of  the  Phase  III  tests  is  shown  in  Table  IV.  Because  of  the  indicated 
difficulties,  quantitative  data  was  acquired  only  at  the  PF3  condition.  The  TVA 
shutdowns,  caused  by  electrical  power  shortages  in  the  AEDC  area  due  to  prevailing 
weather  conditions,  were  particularly  serious.  Since  the  tunnel  shift  started  at 
9:00  a.m.,  and  approximately  3  hours  were  required  to  bring  the  tunnel  to  test 
condition,  equilibrate  the  model,  and  prepare  the  hot-wire  Cor  data  taking,  the 
early  afternoon  shutdowns  precluded  the  acquisition  of  any  quantitative  data  dur¬ 
ing  the  last  three  days  of  this  entry. 

The  difficulties  encountered  with  the  probe  actuator  and  its  vacuum  enclosure 
were  not  unexpected  since  this  device  had  not  beer  operated  with  air-on  prior  to 
the  Phase  III  test.  Although  the  problems  associated  with  the  actuator  were  even¬ 
tually  corrected,  the  lost  time  slowed  the  progress  of  the  tests.  It  should  be 
pointed  out,  however,  that  the  experience  acquired  with  the  actuator  during  this 
entry  led  to  essentially  trouble-free  operation  during  the  Phase  IV  entry. 

The  data  collected  at  the  PF3  condition  was  of  relatively  poor  quality.  In  parti¬ 
cular,  the  hot-wire  signal  was  obscured  by  a  large  noise  signal  which  would  appear 
suddenly  as  high  amplitude  "bursts"  or  spikes.  It  was  originally  suggested  that 
the  noisy  signal  was  caused  by  vibration  of  the  hot-wire  probe  or  by  variations 
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in  the  flow  associated  with  streamwise  shifting  of  the  transition  region.  During 
the  Phase  IV  entry  the  noisy  signal  was  traced  to  a  defective  or  dirty  switch  and 
a  short  length  of  unshielded  cable.  Measurements  at  the  PF3  condition  were  re¬ 
peated,  therefore,  during  the  Phase  IV  tests.  However,  in  order  to  obtain  an 
estimate  of  the  fluctuation  intensities,  the  FF3  uncorrected  wideband  data  at  the 
x  o'  29  inch  station,  where  the  signal  was  relatively  well-behaved,  was  subjected 
to  a  modal  analysis  which  indicated  peak  temperature  fluctuations  of  approximately 
30%.  During  the  analysis  it  was  found  that  the  hot-wire  was  operated  over  a  very 
limited  range  of  overheat  currents  which  Introduced  large  uncertainties  in  the 
results.  Nevertheless,  the  result  was  consistent  with  those  of  KiBtler  (Reference 

9)  extrapolated  to  the  present  test  condition,  and  of  Fischer,  et  al.  (Reference 

10)  obtained  in  adiabatic  helium  flows  at  Mach  20,  and  later  justified  by  the 
measurements  made  during  the  Phase  IV  entry. 

In  addition  to  the  above  data,  a  considerable  amount  of  information  was  obtained 
in  the  form  of  shadowgraph  photographs  and  qualitative  hot-wire  traverses.  This 
latter  information,  which  was  extremely  helpful  in  planning  the  Phase  IV  entry, 
was  used  to: 

(1)  Establish  the  relative  nature  of  the  boundary- layer  flow  as  a  func¬ 
tion  of  Reynolds  number,  blowing  rate,  and  x-station. 

(2)  Delineate  regions  of  transition  flow. 

(3)  Define  the  boundary- layer  thickness  and  its  x-dependence  as  a 
function  of  test  conditions. 

(4)  Indicate  the  relative  effect  of  blowing  on  the  boundary- layer  flow. 

(5)  Provide  an  indication  of  turbulent  activity  across  the  boundary- 
layer. 

Shadowgraphs  of  the  flow  over  the  cone  at  a  free-stream  Reynolds  number  of  1.7 
x  10^/ft  are  shown  in  Figure  22  for  blowing  rates  of  X  ■  0  and  0.0015.  The  left- 
hand  window  shows  the  forward  portion  of  the  cone  extending  over  the  region  from 
x  ■  4  inches  to  x  ■  21  inches  while  the  right-hand  window  indicates  the  region 
from  x  ■  29  to  x  ■  41  inches.  The  grid  spacing  in  Figure  22  is  4^  inches  with 
the  left-hand  vertical  grid  line  located  at  x  ■  29.2  inches  and  the  right-hand 
vertical  grid  line  at  x  *  38,2  inches.  With  \  ■  0  a  lRminar-like  boundary  flow 
extends  to  an  x-station  in  excess  of  30  inches.  With  blowing  the  boundary  thickens 
by  a  factor  of  2  or  3  with  the  transition  zone  shifting  forward  6-8  incheB.  In 
both  cases  a  high-frequency  wave  structure  which  was  also  detected  by  the  hot¬ 
wire  is  apparent  at  the  outer  edge  of  the  layer  in  the  region  of  laminar  flow.  A 
Bimilar  structure  has  been  observed  by  several  other  investigators  (References 
11  and  12),  although  neither  its  origin  nor  its  effect  in  transition  has  been 
identified.  The  scale  of  the  wavefront  is  about  twice  the  boundary-layer  thick¬ 
ness  so  that  the  corresponding  frequency  can  be  given  approximately  by  u/25, 
which  is  in  agreement  with  other  findings,  The  present  hot-wire  measurements 
reveal  the  additional  fact  that  the  wavefront  gives  rise  to  strong  localized 
fluctuations  which  are  similar  in  magnitude  to  the  peak  fluctuation  intensities 
observed  in  the  turbulent  portion  of  the  boundary- layer . 
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In  order  to  provide  a  complete  picture  of  the  relative  nature  of  the  boundary- 
layer  flow  ac  a  function  of  test  condition,  shadowgraphs  obtained  during  the 
Phase  IV  entry  are  included  here.  Shadowgraphs  obtained  ac  a  free- stream  Rey¬ 
nolds  number  of  2.6  x  10°/ft  (PQ  ■  600  psia)  and  shown  in  Figure  23  for  \  ■  0 
and  X  ■  0.0013,  while  those  obtained  at  a  Reynolds  number  of  3.7  x  106/ft  (PQ 
■  850  psia)  for  the  same  blowing  rates  are  shown  in  Figure  24.  These  photo¬ 
graphs  reveal  similar  features  as  those  observed  at  the  low  Reynolds  number 
condition  and  indicate,  in  addition,  a  thinning  of  the  boundary- layer  and  a 
forward  shift  of  the  transition  region  as  the  Reynolds  number  is  increased. 
Furthermore  the  high-frequency  wave  structure  occurring  et  the  edge  of  the 
laminar  boundary- layer  peraists  over  the  entire  Reynolds  number  range. 

The  relative  effect  of  the  free-stream  Reynolds  number  and  blowing  and  the  boun¬ 
dary-layer  turbulence  le  Indicated  in  Figure  23  which  shows  continuous  traces  of 
the  rms  hot-wire  output  obtained  at  x  ■  29  inches  for  a  constant  moderate  over¬ 
heat  current  and  several  test  conditions.  At  a  Reynolds  number  of  1.7  x  10°/ft 
the  boundary- layer  flow  at  x  ■  29  inches  is  laminar  and  the  hot-wire  signal  is 
localized  at  the  outer  edge  of  the  layer  in  the  vicinity  of  the  high-frequency 
wave  structure  observed  in  the  shadowgraphs.  There  is  no  apparent  turbulent 
activity  within  the  layer  where  the  hot-wire  signal  drops  to  its  frea-stresm 
value.  At  a  blowing  rate  of  X  ■  0.0013,  the  boundary  layer  remains  laminar  but 
is  appreciably  thickened  and,  according  to  pitot  pressure  measurements,  moved 
physically  away  from  the  cone.  Again  the  peak  hot-wire  signal  is  localized  at 
the  edge  of  the  boundary- layer  and  there  is  no  apparent  "turbulence11  within  the 
layer,  particularly  near  the  wall.  When  the  free-stream  Reynolds  number  is  in¬ 
creased  to  3.7  x  106/ft,  the  boundary- layer ,  with  X  ■  0.0015,  is  now  fully  turbu¬ 
lent.  The  peak  hot-wire  signal  occurs  at  approximately  the  mid-point  of  the 
boundary- layer  and  fluctuations  persist  to  the  cone  surface  where  the  raw  wire 
output  is  still  well  above  its  free-stream  level. 

3.4.2  PHASE  IV  TESTS 

The  Phase  IV  Philco-Ford  entry  in  Tunnel  B  was  successfully  completed  with  quan¬ 
titative  boundary- layer  turbulence  surveys  made  at  all  of  the  conditions  listed 
in  Table  III  except  the  x  ■  55  inch  station  for  PF6,  A  50  micro-inch  diameter 
hot-wire  was  used  for  the  measurements  at  350  psia  with  x  stations  of  19,  29  and 
39  inches.  All  other  turbulence  measurements  were  made  using  a  20  uinch  dia¬ 
meter  wire.  After  the  difficulties  causing  the  "noisy"  hot-wire  signals  dea- 
cribed  in  Section  3.4.1  were  corrected,  the  noise  was  reduced  to  a  low  level  and 
signal  quality  was  significantly  improved. 

Although  it  was  originally  intended  to  carry  out  the  initial  tests  of  this  entry 
at  850  psia,  excessive  hot-wire  breakage  was  encountered  at  this  condition.  Con¬ 
sequently,  the  maximum  tunnel  pressure  was  restricted  to  600  psia  which,  never¬ 
theless,  should  provide  data  at  the  rearward  x  stations  that  is  representative 
of  highly  developed  turbulent  boundary- layers .  It  is  now  believed  that  the  un¬ 
usually  high  failure  rate  can  be  attributed  to  the  semi-annual  two-week  down-time 
which  preceded  the  test  entry.  During  this  period  the  nozzle  was  recoated  and 
the  interior  of  the  test  section  was  cleaned  and  repainted.  Since  the  tunnel 
was  not  operated  prior  to  the  test  entry,  there  was  no  opportunity  to  eliminate 
the  particulate  material  which  accumulated  during  the  down-time  or  was  eroded 
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FIGURE  23.  SHADOWGRAPH  OF  G 
AND  P  =  600  psi 
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FIGURE  25.  TYPICAL  VARIATION  OF  WIDEBAND  HOT-WIRE  SIGNAL  AT  FIXED 
OVERHEAT  CURRENT  ACROSS  THE  BOUNDARY  LAYER  SHOWING 
EFFECT  OF  REYNOLDS  NUMBER  AND  BLOWING  RATE. 
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from  the  newly  treated  surfaces.  It  Is  recommended,  therefore,  that  subsequent 
hot-wire  tests  should  not  be  scheduled  to  immediately  follow  a  tunnel  down-time. 

During  these  tests  it  was  also  observed  that  wires  with  a  large  aspect  ratio 
(i.e.,  1/ d>  200)  have  poor  survival  times  in  the  flow,  This  accounts  for  the 
lack  of  turbulence  data  with  blowing  at  the  x  ■  55  inch  station.  These  tests 
were  performed  at  the  end  of  the  entry  when  the  short  aspect  ratio  wires  were  no 
longer  available.  Subsequent  hot-wire  tests  should  be  carried  out,  therefore, 
using  wires  with  aspect  ratios  in  the  range  from  1.50  to  180.  The  present  test 
results  confirm  that  in  spite  of  the  large  end  loss  effectB  associated  with  short 
aspect  ratios,  the  unit  Reynolds  number  in  Tunnel  B  is  sufficiently  large  that 
such  wires  have  an  adequate  sensitivity  for  turbulence  measurements. 

It  should  be  pointed  out  that  the  probe  retraction  and  traversing  mechanism  oper¬ 
ated  in  a  satisfactory  manner  and  was  effectively  utilized  in  the  installation 
of  the  fifteen  hot-wiras  used  during  the  experiment.  Most  important,  the  interest 
and  cooperation  of  the  ARO  personnel  assigned  to  this  task  was  instrumental  in  the 
success  of  the  test  program. 

An  example  of  the  spectral  content  of  the  turbulent  fluctuations,  as  well  as  an 
indication  of  the  exceptional  quality  of  the  signals,  is  presented  in  Figure  26 
which  was  obtained  by  playback  of  the  hot-wire  output  recorded  at  x  “1.9  and  29 
inches  for  the  PF5  condition  through  the  wave  analyzer  and  to  an  x-y  plotter. 

This  figure  shows  the  frequency  resolved  signals  corresponding  to  the  noise  and 
the  ten  overheat  currents  for  selected  positions  within  the  boundary- layer .  At 
x  *  19  inches,  the  maximum  signal  occurs  at  y/6  r*  0.6  where  it  is  dominated  by  a 
175  kc  component  which  corresponds  closely  to  the  160  kc  frequency  characteristic 
of  the  wave  structure  observed  near  the  edge  of  the  laminar-like  boundary-layer. 

The  hot-wire  signal  quickly  diminishes  to  near  zero  as  both  the  wall  and  the  free- 
stream  ara  approached.  In  contrast,  at  x  ■  29  inches,  the  boundary-layer  is  fully 
turbulent  and  relatively  Large  fluctuations  occur  throughout  the  entire  boundary- 
layer  with  a  peak  occurring  at  y/6  0,25. 

As  an  assist  in  validating  the  recorded  data,  the  AEDC  computer  was  programmed  to 
carry  out  a  wideband  modal  analysis  using  erm8  values  uncorrected  for  the  frequency 
response  of  the  electronics.  The  method  of  analysis  was  baaed  on  the  HOTWIRE  pro¬ 
grams  described  in  Reference  1  which  were  developed  for  reduction  of  similar  houn¬ 
dary-layer  data  acquired  in  the  JPL/HWT.  This  on-site  capability  permitted  the 
calculation  of  the  y  distribution  of  modally  resolved  fluctuation  intensities  for 
boundary- layer  surveys  performed  at  five  x  stations  prior  to  the  end  of  the  Phase 
IV  entry.  In  addition,  results  for  three  more  surveys  were  available  within  sev¬ 
eral  days  following  the  termination  of  the  test.  The  results  of  the  modal  ana¬ 
lysis,  which  at  this  stage  must  be  considered  a  preliminary  estimate,  yielded 
fluctuation  intensities  consistent  with  expectations  based  on  other  experiments, 
verified  that  the  hot-wires  were  operated  with  an  adequate  overheat  range,  and  In¬ 
dicated  in  most  cases  the  exceptional  quality  of  the  raw  signals.  Thu  latter  is 
illustrated  in  Figure  27,  where  mode  diagrams  for  eight  positions  in  the  boundary- 
layer  at  the  x  “  39  inch  station  for  the  PF3  condition  have  been  plotted.  There 
is  relatively  little  scatter  in  the  data,  represented  by  the  open  circles,  and 
the  least  squares  fit  to  the  theoretical  hyperbola,  Indicated  i>y  the  solid  lines, 
is  extremely  good. 
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y/6  =  0.10 


FIGURE  26.  TYPICAL  SPECTRA  OF  HOT-HIRE  SISALS  FOR:  (A)  TBRBI 
BOUNDARY  LAYER,  AND  (B)  "LAMINAR  LIKE"  BOOBDARY  U 
ffAm  CCBVE  CORRESPONDS  TO  FIXED  OVERHEAT  C WREST. 


;  DIAGRAMS  OF  HOT-WIRE  DATA  FOR  PQ  =  350  psi 

STATION  =  39  INCHES.  OPEN  CIRCLES  INDICATE 
AND  SOLID  LINES  REPRESENT  CURVE  FITS  USED 
.JACTITATION  INTENSITIES. 


The  y-diatributiona  of  temperature  and  velocity  fluctuation  Intensities  (for  the 
purpose  of  data  validation,  the  pressure  fluctuations  were  assumed  to  be  negligible) 
obtained  from  the  on-site  analysis  are  plotted  in  Figures  28  and  29  for  P0  ■  350 
and  600  psia,  respectively.  For  the  low  tunnel  pressure  in  particular,  these  fig¬ 
ures  depict  the  development  of  turbulence  in  the  boundary- layer  growing  over  a  5° 
cone  at  Mach  8  with  an  equilibrium  wall  condition.  They  indicate  further  that 
blowing  at  X  ■  0.0015  increases  the  peak  fluctuations  only  slightly  over  the  non- 
blowing  case,  and  more  significantly,  that  the  temperature  fluctuations  in  both 
cases  are  as  large  as  30-40%.  Clearly,  fluctuations  of  this  magnitude  cannot  be 
handled  by  the  linear  theory  which  forms  the  basis  of  the  conventional  mode  ana¬ 
lysis  and  final  reduction  of  the  data  will  require  a  reformulation  of  the  basic 
treatment  of  the  problem. 

3.5  STATUS  AND  PLANS 

Final  calculation  of  mean  flow  profiles  by  AEDC  using  mean  flow  measurements  made 
by  Philco-Ford  are  in  progress.  The  preliminary  modal  analysis  carried  out  on¬ 
site  are  being  checkad  in  order  to  delineate  potential  problem  areas  during  the 
final  analysis.  Calculation  of  the  hot-wire  sensitivity  coefficients  needed  for 
the  computer  program  DRA2E,  described  in  Reference  1,  and  digitization  of  the  hot¬ 
wire  fluctuation  signal!  which  is  required  for  frequency  resolved  analysis  of  the 
data,  will  be  Initiated  shortly.  In  view  of  the  large  fluctuation  intensities 
discovered  so  far,  the  hot-wire  theory  will  be  re-examined  and  methods  to  account 
for  non-linear  disturbances  will  be  sought. 
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FIGURE  29.  RESULTS  OF  PRELIMINARY  ANALYSIS  SHOWING  Y 
TEMPERATURE  FLUCTUATIONS  AT  P„  =  600  psi. 


SECTION  IV 


CORRELATION  OF  BOUNDARY- LAYER  FLUCTUATION  RESULTS 


In  various  tasks  under  the  present  contract  results  have  been  obtained  on  the 
distribution  of  turbulence  in  supersonic  and  hypersonic  boundary-layers,  Al¬ 
though  the  test  design  for  each  task  aims  at  simulating  specific  flight  condi¬ 
tions  the  simulation  is  always  incomplete  because  of  state-of-the-art  limita¬ 
tions.  In  this  section  we  shall  attempt  to  see  if  a  coherent  picture  emerges 
from  comparison  of  the  available  results  with  each  other. 

4.1  AVAILABLE  EXPERIMENTAL  DATA 

The  available  experimental  data  are  shown  on  Table  V.  Criteria  for  selecting 
data  from  what  is  available  in  the  literature  were: 

(1)  The  test  conditions  must  be  completely  controlled. 

(2)  Sufficient  data  points  must  exist  to  present  a  coherent  picture, 

(3)  The  measuring  technique  must  be  reliable. 

The  data  shown  cover  the  range  from  Mach  0  to  nearly  10;  the  low-speed  data  are 
included  for  purposes  of  comparison.*  The  momentum  Reynolds  number  Rq  ia  also 
spread  widely  (from  103  to  nearly  4  x  10/4).  Most  of  the  data  are  adiabatic  (or 
nearly  so)  with  no  mass  injection.  The  TBL-C  results  shown  are  still  tentative 
since  the  data  have  not  yet  been  fully  reduced;  they  are  taken  from  Section  III 
of  this  report.  Also,  the  WEM  and  TBL-A  data  are  taken  from  Sections  II  and  VI 
of  this  report. 

4.2  VELOCITY- FLUCTUATION  CORRELATIONS 

Since  there  is  no  theory  available  the  fluctuation  distributions  as  a  function 
of  y/6  can  be  presented  empirically  correlated  into  some  sort  of  a  "curve "  versus , 
say,  Mach  number;  or  it  can  be  further  fortified  into  a  non-dimensional  grouping 
which  "collapses  the  data  into  one  curve"  including  other  parameters  such  as  heat 
transfer . 

Figure  30  shows  the  simplest  non-dimensional  correlation  u'/u  vs,  y/6  where  u  ,1s 
the  local  average  flow  speed.  There  generally,  a  tendency  to  decrease  u'/u 
as  M  ia  ^increased.  When  the  velocity  ue  at  the  boundary- layer  edge  is  used  in¬ 
stead  of  u  the  decrease  in  the  fluctuation  intensity  is  not  obvious  (for  the 


*  There  appear  to  be  several  sets  of  low-speed  fluctuation  data  which  fulfill 
the  above  criteria.  The  data  of  Townsend  (Reference  1'j)  and  KlebRnoff  (Ref¬ 
erence  14)  are  typical  of  these. 
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typical  position  y/6  *  0.5)  until  hypersonic  Mach  numbers  are  reached,  aa  shown 
on  the  right  in  Figure  30.  Nearer  the  wall,  however,  the  correlation  u'/ue(y/6) 
shows  a  continuous  decrease  of  the  fluctuation  magnitude  as  the  Mach  number  de¬ 
parts  from  zero. 

In  his  paper  (Reference  9)  Kistler  follows  the  practice  of  the  low-speed  experi¬ 
menters  to  normalize  u1  with  the  "shear  velocity" 


The  u'/u  correlation,  shown  on  Figure  31,  still  shows  a  strong  Mach  number  de- 
pendence^  even  accounting  for  differences  in  Rg  or  the  temperature  ratio  w 
It  remains  to  be  seen  if  the  correlation  improves  (i.e.,  if  the  curves  collapse) 
when  the  proper  distance  parameter 
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0 


dy 


is  considered  in  lieu  of  the  physical  distance  y/6  shown  on  Figure  31.  The  rea¬ 
soning  for  using  the  u'/uT  vs.  y+  correlation  is  its  well-known  validity  in  col¬ 
lapsing  mean  flow  data  for  high-speed  boundary- layers  (Reference  2), 

Another  possible  correlation  comes  from  Morkovin's  hypothesis  (The  "Dynamic  Equi¬ 
librium  Hypothesis",  Reference  15)  that  the  "universal"  form  of  the  fluctuations 
is  restricted  to  the  turbulent  energy  pu as  opposed  to  the  turbulence  velocity 
u'  alone.  The  proper  parameter  to  plot,  then,  is  the  quantity  (u'/uT) .(p/pw)%. 
This  has  been  attempted  in  Figure  32,  for  which  some  data  (Reference  16)  other 
than  those  discussed  so  far  are  also  shown.  There  is  indeed  some  improvement  over 
the  method  of  plotting  of  Figure  31  (where  the  density  factor  s*'  had  been 

omitted)  in  the  outer  half  of  the  boundary- layer .  Near  the  wa..  ..ever,  this 
latest  plotting  method  still  shows  scatter  by  a  factor  of  more  th.  2  In  the  data, 
Also,  the  turbulent  energy  is  a  more  complex  expression  than  juBt  pu‘ 2  which  is, 
in  fact,  probably  smaller  numerically  than  the  omitted  terms.  It  is  thought 
necessary  to  look  into  this  matter  further,  i.e.,  to  improve  on  Morkovin's  cor¬ 
relation  by  plotting  the  complete  energy  term  vs.  y+. 


4.3  TEMPERATURE-FLUCTUATION  CORRELATIONS 


The  local  wideband  temperature  f luctuations_have  been  plotted  on  Figure  33  nor¬ 
malized  with  the  local  average  temperature  T.  The  graph  on  the  left  docs  not 
Include  the  tentative  TBL-C  data  (see  Section  III),  but  these  results  are  in¬ 
cluded  in  the  cross-plot  (at  y/6  **  0.5)  shown  on  the  right.  It  is  seen  that  the 
adiabatic  or  near-adiabatic  tests  show  a  very  rapjld  increase  of  the  fluctuation 
intensities  with  Mach  number.  Beginning  from  T'/T  ■  0  at  M  ■  0  (by  definition) 
these  fluctuations  have  increased  to  about  22%  at  M  »  7.  It  must  also  be  noted 
that  the  qualitative  data  of  Wagner  and  his  co-workers  (Reference  17),  excluded 
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FIGURE  31.  LONGITUDINAL  WIDEBAND  TURBULENCE 

LAYERS  (uT  CORRELATIONS ) . 


FIGURE  32.  CORRELATION  OF  WIDEBAND  LONGITUDINAL  VELOCITY  FLUCTUATIONS 
BASED  ON  FRICTION-VELOCITY  SCHEME  FOR  COMPRESSIBLE  FLOW. 


here  because  of  the  criteria  of  Section  4.1,  have  shown  T'/T  of  order  100%  at 
M  *  20.  This  result  of  Wagner's  had  been  earlier  subject  to  suspicion  because 
it  was  obtained  with  a  linearized  technique  which  is  invalid  for  such  large 
fluctuations.  Nevertheless,  in  view  of  the  results  of  Figure  33,  it  can  at  least 
be  stated  that  at  M  ■  20  the  T'/T  must  be  very  large  indeed.  On  the  other  hand, 
the  "cooled  wall"  point  of  Figure  33  (from  TBL-A,  Reference  2  and  Section  II  of 
this  report)  shows  a  very  strong  effect  of  hast  transfer  on  T'/T. 

An  attempt  to  better  correlate  the  results  of  Figure  33  was  made  by  modifying 
slightly  Kistler's  original  suggestion  (Reference  9).  This  suggestion  is  that 
the  fluctuations  of  a  passive  scalar  such  as  temperature  are  "driven"  by  the  tem¬ 
perature  difference  T0  -  Te,  i.e.,  the  total  minus  the  static  stream  temperature. 
Actually,  the  driving  force  is  the  temperature  drop  Tr  •  T  across  the  boundary- 
layer,  In  adiabatic  flows  T-  is  very  nearly  T0,  but  with  neat  transfer  Tr  is  the 
wall  temperature  Tw.  Also,  T  is  itself  normalized  with  the  "average"  temperature 
(Tr  +  Tc)/2,  so  that  finally  the  resulting  quantity  plotted  is 
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versus  y/6, 


where  Tv  is  actually  the  wall  temperature.  The  results,  lesB  the  TBL-C  data,  are 
plotted  on  Figure  34.  The  scatter  from  the  data  of  one  experiment  to  the  other 
is  considerably  reduced  in  this  manner.  The  TBL-A  results  are  plotted  in  two  dif¬ 
ferent  ways,  ona  by  assuming  Tr  ■  T0  and  the  other  Tr  ■  Tw.  The  latter  gives  bet¬ 
ter  agreement,  implying  that  the  correlation  is  beBJ  when  the  wall  temperature  is 
indeed  chosen.  The  importance  of  this  finding  1b  that  cases  with  heat  transfer 
can  be  taken  care  of  in  this  way. 

4.4  CONCLUSIONS  FROM  THE  CORRELATION  STUDY 

On  the  basis  of  the  experimental  data  studied  so  far  the  following  conclusions 
are  drawn : 

(a)  There  is  a  tendency  of  increasing  Mach  number  to  decrease  the  longi¬ 
tudinal  velocity  fluctuations. 

(b)  The  best  correlation  is  obtained  by  using  the  "dynamic  equilibrium 
hypothesis"  of  plotting  (u  ' /u<r) •  (p/pw)^  .  However  the  correlation  1b 
good  only  in  the  outer  half  of  the  layer;  it  is  invalid  near  the 
wall  if  the  physical  distance  y/6  is  used, 

(c)  The  static  temperature  fluctuations  T'/T  appear  better  correlated  if 
plotted  versus  2(TW  -  Te)  f  (y/6)/(Tw  +  Te) .  This  scaling  formula  pre¬ 
dicts  that0T'/T  increase  sharply  with  M  and  decrease  greatly  when  the 
wall  is  coaled, 

(d)  Provided  that  the  Layer  thickness  is  much  smaller  than  the  trans¬ 
verse  curvature  there  is  no  effect  of  geometry  (flat  plate  versus 
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FIGURE  34.  CORRELATION  OF  WIDEBAND  TEMPERATURE  FLUCTUATIONS 
SUPERSONIC/HYPERSONIC  BOUNDARY  LAYERS. 


cone),  and  there  seems  to  be  no  great  effect  of  Reynolds  number. 
The  most  Important  influence  found  is  that  of  the  heat  transfer 
on  the  T'/f. 


4.5  STATUS  AND  PLANS 

From  the  numerical  data  presented  and  from  the  above  conclusions  it  i»  now  pos¬ 
sible  to  estimate  the  fluctuation  Intensity  u'/u  and  T'/T  for  boundary- layers  in 
the  range  0  <  M  <  10.  However  the  conclusions  are  tentative  because  the  data 
base  is  very  small',  each  test  discussed  involves  changes  in  at  least  two  para¬ 
meters  compared  with  any  other  test.  Therefore  it  is  not  certain  which  of  the 
two  flow  parameters  create  the  changes  noted  between  any  two  experiments.  Some 
parameters  (e.g.,  heat  transfer)  are  present  in  only  one  experiment,  while  some, 
such  as  mass  addition  or  pressure  gradient,  have  not  as  yet  been  covered. 

In  the  meantime,  it  is  planned  to  continue  the  effort  to  correlate  the  data  along 
the  following  lines: 

(1)  Forthcoming  data  from  the  current  Philco-Ford  experiments  will  be 
analyzed  and  plotted  in  a  manner  such  as  that  of  Figures  30  through 
34.  These  data  Include: 

(a)  TBL-C  data  which  include  wider  variations  in  Rg  (such  as  ob¬ 
tained  by  using  the  profiles  at  different  x-distances  from 
the  cone  tip)  and  mass  addition  to  the  extent  m/peue  ■  0.0015, 

(b)  WEM  data  which  also  include  wider  variations  of  Rg.  The  pro¬ 
file  at  the  cone  shoulder  at  condition  C  (P0  -  530  mm  Hg,  Bee 
Section  VI)  will  be  available  shortly. 

(2)  The  velocity  fluctuation  data  of  Figure  32  will  be  re-plotted  to 
change  the  abscissa  to  y+,  This  will  involve  modifying  other 
worker's  data  into  this  coordinate  system. 

(3)  The  velocity  fluctuations  will  be  further  plotted  with  coordinates 
which  include  the  total  turbulent  kinetic  energy. 
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SECTION  V 


MEASUREMENT  OF  SHEAR  STRESSES  IN  A  HYPERSONIC  BOUNDARY- LAYER 


5.1  PURPOSE 

The  purpose  of  this  task  1b  the  measurement  of  the  turbulent  shear  stress  (the 
so-called  "Reynolds  stress")  in  a  hypersonic  boundary-layer. 

5.2  SUMMARY  OF  PREVIOUS  WORK 

Work  on  this  problem  was  initiated  in  the  fall  of  1971  with  the  development  of 
techniques  for  fabricating  and  calibrating  the  X-array  hot-wire  anemometers  used 
for  the  shear  stress  measurements  and  was  temporarily  halted  in  April  1972  fol¬ 
lowing  the  completion  of  the  experimental  program.  Tests  were  performed  in  the 
JPL  Hypersonic  Wind  Tunnel  under  the  same  conditions  used  in  the  TBL-A  experiment 
(References  1  and  2).  This  provided  the  advantage  of  making  shear  stress  measure¬ 
ments  in  a  boundary- layer  whose  mean  and  fluctuating  flow  had  already  been  well- 
defined,  A  total  of  eight  10,000  foot  tapes  of  recorded  data  were  acquired  repre¬ 
senting  about  60  y  positions  in  the  boundary- layer .  Of  these  data  roughly  one-third 
Is  suspected  to  be  of  inferior  quality  due  to  excessive  scatter  in  the  measured  wire 
resistance,  leaving  therefore  nearly  40  points  of  useful  information. 

In  addition  to  the  measurements,  equations  governing  the  "general"  behavior  of  X- 
probes  in  compressible,  turbulent  flows  were  derived  (assuming  negligible  pressure 
fluctuations)  providing  the  means  for  deducing  the  Bhear  stress  terms  from  the 
hot-wire  signals.  Preliminary  calculations  indicated  that  the  fluctuations  u7^, 

T7"^  and  u  rT7  are  numerically  similar  to  those  found  with  a  single  wire  in  the  TBL-A 
experiment;  thereby  establishing  confidence  in  the  validity  of  the  remaining  un¬ 
knowns,  Trv* ,  u  V  and  v^,  In  addition,  approximate  calculations  of  the  latter 
parameters  were  carried  out  at  three  positions  on  the  boundary- layer .  The  v1^  was 
found  to  lie  numerically  below  u7^  as  expected,  while  the  correlations  u 'v 1  and 
v'T'  were  found  to  be  negative  as  anticipated  near  the  wall  but  to  change  sign 
near  the  outer  edge.  It  is  not  known  whether  the  latter  result  is  real  or  a  con¬ 
sequence  of  the  incompleteness  of  the  preliminary  analysis. 

5.3  PRESENT  STATUS 

Since  the  TBL-A  experiments  indicated  the  importance  of  pressure  fluctuations  at 
hypersonic  speeds,  the  X-probe  theory  developed  previously  hns  been  extended  to 
imlude  this  effect.  Details  of  this  formulation  are  presented  in  Appendix  B, 

In  the  case  of  negligible  pressure  fluctuations,  the  X-probe  theory  leads  to  six 
linear  simultaneous  equations  of  sufficiently  simple  form  to  be  solved  by  hand 
if  necessary.  With  finite  pressure  fluctuations,  several  restrictive  assumptions 
concerning  the  correlation  between  the  sound  and  other  modes  are  required  and 
while  six  governing  equations  are  again  obtained,  in  this  case  the  equations  are 
non-linear  and  cannot  be  solved  without  the  aid  of  a  computer.  For  this  reason, 
and  since  the  results  of  TBL-A  Indicated  that:  the  magnitude  of  velocity  and  tem¬ 
perature  fluctuations  obtained  from  a  no-sound  assumption  are  Increased  only 
slightly  above  their  corresponding  values  when  relatively  large  pressure 
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fluctuations  are  accounted  for,  it  is  planned  to  reduce  the  X-probe  data  in  two 
stages.  First,  a  solution  will  be  found  by  assuming  a  negligible  sound  field, 
thus  providing  an  opportunity  to  compare  u'"?,  0^  and  u  '0 '  with  the  no-sound  re¬ 
sults  of  TBL-A  and  to  assess  the  relative  quality  of  the  X-probe  data.  When  this 
is  accomplished  the  more  complex  solution  allowing  finite  pressure  fluctuations 
will  be  carried  out. 

Calculation  of  the  hot-wire  sensitivity  coefficients,  which  are  needed  to  inter¬ 
pret  the  X-probe  data,  has  also  been  initiated.  This  calculation  requires  a 
knowledge  of  the  heat  transfer  characteristics  of  the  wire  which  is  usually  ac¬ 
quired  from  a  flow  calibration  yielding  the  Nusselt  number  -  Reynolds  number  and 
recovery  factor  -  Reynolds  number  relations  of  the  probe.  The  flow  calibration 
is  generally  accomplished  by  locating  the  wire  in  the  free  stream  of  the  wind- 
tunnel,  and  prior  to  these  tests,  the  calibration  data  was  obtained  by  hand-record¬ 
ing  several  1,  e  pairs  in  order  to  find  the  zero  current  R  and  dR/di2  values  used 
to  evaluats  Nu  and  T].  From  the  overheat  traverses  carried  out  during  the  boundary- 
layer  survey,  the  local  Nu  can  be  found  and  used  with  the  heat  transfer  calibra¬ 
tions  to  determine  the  corresponding  derivatives  d  log  Nu/d  log  Re  and  d  log  7]/d  log  Re 
appearing  in  the  expression  of  the  sensitivity  coefficients  (see  for  example  Refer¬ 
ence  17). 

Because  of  drift  in  tunnel  conditions  during  the  time  required  for  calibration  and 
scatter  in  the  hand-recorded  i,  e  data,  significant  uncertainty  la  frequently  in¬ 
troduced  in  the  heat  transfer  calibration.  In  order  to  improve  the  calibration, 
an  analog  circuit  was  designed  and  first  used  during  the  X-probe  tests.  The  cir¬ 
cuit  receives  as  input  instantaneous  i,  e  values  and  outputs,  directly  to  an  X-Y 
plotter,  the  corresponding  R  and  i2 .  The  calibration,  at  a  given  runnel  condition, 
is  obtained  in  a  matter  of  seconds  by  increasing  the  current  from  zero  to  a  pre¬ 
selected  maximum  value  by  means  of  a  single  potentiometer  acting  as  a  series  re- 
sictor  in  the  hot-wire  circuit. 

To  check  the  performance  of  the  analog  circuit  during  its  first  usage  with  X-probe 
No.  4,  calibration  data  was  also  obtained  by  hand-recording  i,  e  pairs.  The  re¬ 
sults  are  shown  in  Figures  35  and  36  where  the  analog  calibration  is  indicated  as 
triangles  and  the  hand-recorded  data  as  open  circles.  The  solid  lines  drawn 
through  the  data  points  represent  a  least-squares  fit  to  the  analog  calibration 
data.  A  alight  but  consistent  difference  is  observed  for  both  wires  except  for 
the  Nu-Re  relation  for  wire  4-2A,  where  there  Is  a  large  divergence  between  the 
two  calibration  procedures.  This  discrepancy  was  resolved  deriving  a  calibration 
using  local  Nusselt  numbers  determined  from  the  boundary- layer  survey  and  local 
Reynolds  numbers  and  total  temperatures  from  the  mean  flow  measurements  of  TBL-A. 

This  "boundary- layer"  calibration  is  depicted  as  squares  in  Figures  35  and  36. 

The  resulting  Nu-Re  relation,  Figure  35,  is  seen  to  be  in  excellent  agreement 
with  the  analog  calibration  for  both  wires  (except  for  the  high  Re  data  for  wire 
4-2A),  However,  the  corresponding  T|-Re  relation,  Figure  36,  shows  poor  agreement 
with  the  free-stream  results,  particularly  for  wire  4-2A  which  exhibits  a  be¬ 
havior  which  cannot  be  rationalized  on  physical  grounds.  The  "boundary- layer" 
calibration  demonstrates  the  superiority  of  the  analog  technique  for  obtaining 
an  accurate  Nusselt-Reynolds  relation  and,  until  an  explanation  Is  evolved  for 
concerning  the  recovery  factor  data  derived  from  the  boundary- layer  measurements, 
the  free-stream  calibration  of  this  parameter  will  be  considered  to  provide  a 
reliable  T|-Re  correlation, 
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FIGURE  35.  NUSSELT  NUMBER  -  REYNOLDS  NUMBER  CALIBRATION  FOR  X- FRO BE 


FIGURE  36.  RECOVERY  FACTOR  -  REYNOLDS  NUMBER  CALIBRATION  FOR 


A  similar  calibration  difficulty  was  encountered  with  X-probe  No.  11,  for  which 
only  the  analog  circuit  was  used.  For  wire  11-1A,  the  current  range  was  insuf¬ 
ficient  to  provide  accurate  values  of  dR/di2  and  the  resulting  Nusselt  calibra¬ 
tion,  shown  for  both  wires  11- 1A  and  11-2A  as  the  triangles  in  Figure  37,  was  sus¬ 
pect.  Again  a  "boundary-layer"  calibration  was  carried  out  which  differed  consid¬ 
erably  from  the  analog  calibration  of  wire  11-1A  although  good  agreement  with  that 
of  wire  11-2A  was  obtained.  This  was  interpreted  again  as  indicating  the  validity 
of  the  analog  calibration  technique  provided  a  sufficient  current  range  is  used. 

The  T|-Re  calibration  for  the  two  wires,  Figure  38,  showB  the  same  behavior  for  the 
boundary- layer  data  exhibited  by  X-probe  No.  4,  although  the  data  for  wire  11-2A 
approaches  much  closer  the  free- stream  analog  results,  This  ImplleB  that  the  free- 
stream  calibration  of  T]  versus  Ro  is,  indeed,  a  reliable  indication  of  the  recovery 
factor  characteristic  of  the  probe. 

As  a  result  of  the  preceding  exercise,  an  alternate  method  of  data  reduction  was 
evolved  to  smooth  the  small  experimental  errors  in  the  Nusselt  number  measurement 
that  are  apparent  in  Figures  35  and  37.  In  the  TBL-A  program,  the  logarithmic 
derivatives  of  Nu  and  T|  were  determined  by  calculating  the  zero  current  R  and  dR/di2 
from  the  local  overheat  traverse,  using  an  iterative  technique  to  find  Re,  Nu,  and 
T) ,  and  evaluating  the  derivatives  from  the  curve-fitted  expressions  to  the  heat 
transfer  calibrations.  In  the  present  case,  however,  it  has  been  decided  to  use 
the  well-defined  mean  flow  profiles  to  specify  the  local  Re,  so  that  the  deriva¬ 
tives  can  be  calculated  directly  from  the  calibration  equations,  thereby  ensuring 
a  smooth  monatomic  increase  (or  decrease)  in  the  value  of  theBe  quantities.  This 
should  lead  to  a  substantial  improvement  in  the  accuracy  of  the  sensitivity  coef¬ 
ficients  and,  eventually,  of  the  modal  analysis. 

The  following  tasks  remain  before  final  analysis  of  the  X-probe  data  can  be  car¬ 
ried  out! 

(1)  Complete  calculation  of  sensitivity  coefficients. 

(2)  Digitize  hot-wire  signals  for  computer  processing. 

(3)  Modify  and  extend  computer  program  DRA2E  to  accommodate  equations 
and  logic  required  for  X-probe  analysis. 
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FIGURE  37.  NUSSELT  NUMBER  -  REYNOLDS  NUMBER  CALIBRATION 


RECOVERY  FACTOR  -  REYNOLDS  NUMBER  CALIBRATION  FOR 
X- PROBE  NO.  4. 


SECTION  VI 


THE  NEAR  WAKE  OF  A  CONE 


6.1  SUMMARY  OF  PREVIOUS  WORK 

The  motivation,  experiment  design  and  teat  results  on  the  mean  flow  for  the  near¬ 
wake  experiments  have  already  been  reported  in  Reference  5  (pp  3-44).  A  summary 
of  these  will  be  given  in  the  present  section. 

The  near-wake  study  requires  the  measurement  of  the  mean  and  turbulent  flow  in 
the  near-wake  of  a  slender  cone.  Since  theory  (Reference  18)  expects  the  major 
effects  to  depend  on  Reynolds,  rather  than  on  Mach  number,  the  experiment  is 
done  in  the  Mach  3.0  flow  of  the  Philco-Ford  Supersonic  Wind-Tunnel,  Accordingly, 
after  examination  of  various  model-suspension  schemes,  a  S-degree  half-angle  cone 
with  a  one-inch  base  diameter  was  suspended  in  the  test  section  by  wires.  A 
"trip"  consisting  of  a  0-ring  provides  early  transition  in  the  model  boundary- layer, 
at  the  same  time  giving  an  acceptably  thin  layer  at  the  shoulder  (6/D~  0.08).  It 
was  found  that  the  tripped  boundary-layer  profiles  were  very  similar  to  those  of 
a  naturally- turbulent  (untripped)  layer,  so  that  the  trip  can  be  considered  equi¬ 
valent  to  a  Reynolds  number  increase. 

Much  effort  was  devoted  to  the  choice  of  number,  size  and  location  of  the  wire 
supports.  A  total  of  seven  wires  ware  finally  used  in  a  way  allowing  their  added 
use  for  pitch  and  yaw  control.  Although  the  presence  of  each  wire  could  be  de¬ 
tected  by  a  sensor  in  its  wake,  the  effect  of  the  supports  on  the  cone  wake  itself 
was  judged  minimal.  This  conclusion  was  arrived  at  by  first  making  measurements 
of  the  cone  base  pressure  P^,  the  location  xn  of  the  "neck"  and  the  wake  thickness 
dy.  The  magnitude  of  Pt,,  xn  and  d„  so  found  agrees  very  well  with  those  of  te3ts 
done  elsewhere,  including  tests  with  magnetically-suspended  bodies. 

Special  pitot  and  static-pressure  probes  were  constructed  and  were  satisfactorily 
checked  for  insensitivity  to  pitch  and  yaw  effects.  A  complete  "matrix"  of  mean- 
flow  data  were  taken  with  these  probes  for  the  conditions  of  Figure  39.  The  data 
were  collected  by  an  automated  system  producing  IBM  cards  in  real  time,  and  were 
reduced  to  give  the  local  mean  velocity,  temperature,  pressure,  density,  total 
temperature  and  Mach  number  at  each  point;  the  relevant  computer  programs  are 
called  WEM-I  and  WEM-II.  A  limited-distribution  report  (Reference  19)  has  been 
issued  which  contains  plots  of  all  the  mean  quantities  for  the  conditions  of 
Figure  39.  Examples  of  these  are  also  shown  in  Figures  16,  21  and  22  of  Reference 
5.  The  major  conclusions  drawn  from  the  work  of  the  preceding  period,  as  described 
Just  above,  were! 

(1)  The  experiment  design  goals  were  met:  the  test  conditions  ranged 
from  the  so-called  "pre-BLT"  condition  (laminar  boundary-layer,  at 
Re  D  ■  69,000)  to  the  "post-BLT"  conditions  (turbulent  layer, 

Re“*n  >  100,000), 

oo  y  w 

(2)  The  initial  conditijns  (i.e.,  the  boundary- layer  profiles)  «t  x  ■ 

0  were  normal. 
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FIGURE  39.  EXPERIMENT  MATRIX  FOR  THE  WEM  (NEAR  WAKE)  TESTS 


(3)  The  wake  velocity  and  temperature  showed  distinct  changes  from 
the  pre-BLT  to  the  post-BLT  conditions. 

(4)  There  were  substantial  pressure  gradients  in  x  <  6D;  the  peak  pres¬ 
sure  is  encountered  beyond  the  neck. 

(5)  The  wake  diameter  Increases  with  x  for  pre-BLt  but  is  constant 
for  post-BLt. 

6.2  FURTHER  REDUCTION  OF  THE  MEAN-FLOW  DATA 

In  the  present  reporting  period  the  mean-flow  data  reported  in  last  period's  re¬ 
port  (Reference  3)  were  further  put  into  a  form  suitable  for  comparison  with 
theories  and  other  experiments  and  for  use  as  generalized  "rules"  in  flight  ap¬ 
plications.  This  was  done  by  computer  program  DRB3E(WEM-III) .  An  outline  of 
these  operations  will  be  given  below. 

To  test  whether  the  wake  approached  gasdynamic  equilibrium  the  wake  radial  co¬ 
ordinate  y  was  transformed  into  the  Howarth  distance 


where  p  is  the  local  density  and  p  the  density  just  external  to  the  wake,  at  a 
distance  ye  from  the  axis;  the  ye  ewas  determined  by  inspection  of  the  radial 
profiles  from  Reference  19.  Using  the  velocity  ue  and  temperature  Te  at  ye,  the 
velocity  and  temperature  "defects" 

u  -  u(y  -  0)  T(y  *  0)  -  T 

w  a  - - - ,  6  s  - - - -  (6.2) 


were  also  computed.  Next,  Che  wake  drag  was  computed  from  the  data  as  follows. 
It  is  known  that  the  drag  integral  at  each  fixed  distance  x  behind  the  cone  is 
given  (Reference  20,  p.  93)  by 


S  n 

wake  drag  ■ 

S1  +S2 

■  viscous  drag  +  pressure  drag 

S1  " 

4h  r  (i 

J  Peue 

-  —  >y  dy 

e 

-  SL(x) 

(6.3) 

S2  " 

-VEf  - 

yM  i  '  L  e 
e 

ijy  yd 

■  $2 (x) 

(6.4) 
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The  integration  limits  for  the  viscous  drag  are,  obviously,  y  ■  0  and  ye. 

The  choice  of  limits  is  not  so  clear  for  the  pressure  drag  S2.  If  the  upper 
limit  extends  beyond  the  cone's  bow  shock  then  S  is  bound  by  the  total  drag  of 
the  cone  rather  than  the  wake  drug.  On  the  other  hand,  S2  cannot  be  completely 
ignored  since  P  j*  Pg,  i.e.,  substantial  pressure  gradients  exist,  as  shown  on 
Figure  21  of  Reference  5.  A  reasonable  compromise  was  chosen  by  computing  the 
pressure  drag  only  within  the  wake  boundary  ye;  therefore,  the  limits  y  ■  0  and 
y  ■  ye  were  also  chosen  for  the  S2  integral.  (Actually,  both  Sj^  and  S2  are 
separately  computed  and  displayed  in  the  program.) 

With  the  defects  w  and  8  and  the  drag  S(x)  known  at  each  x,  the  wake  character¬ 
istic  length  L  (the  "traverse"  scale)  was  computed: 


L(x) 


r  ' 

L  4ttW  . 


(6.5) 


The  similarity  parameter  T|  sa  'y'/L  (i.e.,  the  non-dimensional,  "compressible" 
radial  distance)  was  also  computed  as  well  as  the  non-dimensional  velocity  and 
temperature: 


ui(x,y)  h 


u  -  u(y  ■  0)’ 
6 


T(x,y)  « 


T 


T(y  -  0)  -  T 


(6.6) 


The  corresponding  non-dimensional  axial  distance  is 


/s 

Finally  the  DRB3E  program  computes  the  usua 
"Turbulent  Reynolds  Number"  Rt), 

"Turbulent  Pvandtl  Number"  a.^ 

as  well  as  the  Reynolds  number  Re 

w 

where  M  and  v  are  the  Mach  number  and  la 

p  p 

the  wake. 


(6.7) 


transport  coefficients: 


,  1/2 
6rt 

3/2  - 

w  X 

(6.8) 

4 

(6.9) 

0.4  M 

2 

fcw 

(uc  - 

u(y  -  U))L 

(6.10) 

nav  kinematic  viscosity  Just  outside 
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In  addition,  calculations  of  the  cone  drag  budget  (see  Appendix  C)  and  of  the 
integral  thicknesses  of  the  boundary- layer  were  made  in  order  to  characterize 
the  wake  for  proper  placement  in  the  evolving  matrix  of  laboratory  and  flight 
data. 

6.3  RESULTS  OF  THE  MEAN  FLOW  MEASUREMENTS 
6.3.1  CHARACTERIZING  PARAMETERS 

The  significance  of  the  present  tests  lies  heavily  on  the  ability  to  demonstrate 
the  regularity  of  the  wake  initial  conditions  and  the  capacity  of  these  conditions 
to  generate  meaningful  turbulent  near-wake  phenomena.  Furthermore,  parameters 
such  as  the  cone  drag,  the  momentum  Reynolds  Number  Rq,  etc.,  must  be  known  in 
order  to  match  them  with  the  observed  turbulent  behavior.  The  following  measure¬ 
ments  are  important  in  this  respect! 

(a)  As  already  mentioned  in  Subsection  6.1  end  demonstrated  in  p.  37  of 


Reference  5, 

the  boundary-layer  profiles 

conform 

to  the  expectations 

for  all  test 

conditions  (e.g 

. ,  both  for 

pre-  and 

post-BLT 

conditions) 

The  boundary 

'layer  momentum 

thicknesses  1 

Rg  were  i 

computed  from  the 

measured  profiles  at  the  cone  base  as  fo 

llowa : 

TUNNEL 

CONE 

EXPERIMENT 

PRESSURE  P0 

REYNOLDS  NO. 

e 

D 

6* 

DESIGNATION 

(mm  Hk  aba .) 

(Ran.®) 

(inch) 

!i 

(inch) 

A 

7  30 

178,000 

0.0064 

1,140 

0.0322 

B 

630 

152,000 

0.006 

910 

0.0325 

C 

530 

127,000 

0.0069 

875 

0.0314 

D 

430 

106,000 

0.0069 

730 

0.0296 

E 

330 

81,000 

0.0045 

365 

0.017 

F 

290 

69,000 

0.0042 

290 

0.015 

In  this  table  we  see  that  the  computed  Rg's  confirm  a  posteriori  the 
transition  threshold  observed:  that  is,  the  A-D  experiments  (PQ  >  330 
mm  Hg  abs.)  were  observed  to  have  a  clearly  turbulent  (post-BLT)  boun¬ 
dary-layer,  which  could  have  been  predicted  also  because  their  Rg's  arc 
greater  than  a  few  hundred.  Note,  also,  that  the  "most  turbulent"  test 
condition  produced  (Experiment  A,  PQ  ■  730  mm  Hg  abs.)  involves  an  Rq 
■  1,140,  which  is  probably  much  smaller  than  the  Kg  encountered  in 
flight.  This  point  is  important,  because  Rq  is  a  measure  of  the 
"maturity"  of  the  boundary- layer .  The  present  thinking  is  that  above 
a  certain  Rq  most  boundary- layer  properties  are  Reynolds-number- 
independent  . 

(c)  Another  useful  "benchmark"  In  the  present  tests  is  that,  beyond  a  cer¬ 
tain  distance  behind  the  cone  the  wake  properties  should  tend  toward  the 
familiar  "far  wake"  behavior: 
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w  *-  6 


(6.11) 


--2/3 


Also,  we  expect  the  radical  variations  u  and  T  (i.e.,  Eqna .  (6.b)  to 
fall  Into  a  Gaussian  shape,  that  is 


y\  =  y/L 


k  “  const. 


(6.120 


with  corresponding  asymptotic  formulas  for  the  turbulence  properties. 
These  relations  are  expected  to  hold  not  only  on  theoretical  grounds 
(Reference  21)  but  also  because  they  were  found  valid  during  the  earlier, 
ARPA-sponsored  studies  at  Philco-Ford  of  the  axisymmetric ,  far-wake  be¬ 
havior  at  Mach  3  (References  22  and  23  for  example).  We  therefore  shall 
be  especially  careful  to  recognize  in  the  presentation  of  results  below 
(Subsection  6,3.3.b)  the  relaxation  process  which  will  turn  the  near¬ 
wake  behavior  into  the  far-wake  behavior.  Naturally,  since  the  present 
results  stop  at  x/D  ■  6,  this  asymptotic  behavior  may  not  be  fully  at¬ 
tained.  Any  semblance  of  a  trend,  however,  will  reinforce  the  conviction 
that  the  present  flow  field  begins  (as  discussed  in  Subsection  b.3.3  re¬ 
ferring  to  the  boundary-layer)  and  ends  with  familiar  flow  profiles. 


6.3.2  EXPECTATIONS  FROM  INTUITION  AND  FINSON'S  THEORY 

Before  the  mean-flow  results  are  presented,  the  expectations  of  the  near-wake  be¬ 
havior,  especially  as  enunciated  by  Finson  (Reference  18)  will  bo  discussed  briefly. 
Such  a  discussion  is  needed  for  comparison  with  the  present  experimental  results  to 
see  if  the  present  basic  ideas  for  near-wake  formation  are  valid.  In  this  section 
Finson's  predictions  will  be  simply  referred  to  as  the  "theory". 

The  theory  is  basically  an  extension  of  intuitive  thinking  about  the  relaxation 
process  connecting  the  boundary- layer  at  the  shoulder  and  the  far  wake.  It  says 
that  the  small-scale  turbulence  shed  by  the  body  boundary- layer  persists  for  some 
dozens,  perhaps  hundreds,  of  (drag)  diameters  behind  the  vehicle  without  deter¬ 
ioration  provided  that  its  "turbulence  Reynolds  Number" 


% 


(6.13) 


is  much  higher  than  about  40  m  the  near  wake.  With  small  scales  dominating,  the 
large-scale  motion  characteristic  of  a  typical  turbulent  wake  cannot  develop. 

Thus  the  familiar  "eddy  viscosity"  is  aLso  smaller,  the  spreading  of  the  wake  Is 
inhibited  and  the  "wake  velocity"  decays  very  slowly.  This  creates  a  "plateau" 
in  the  w  versus  x  plot,  whereby  the  wake  velocity  is  nearly  constant  lor  as  long 
as  the  small  scales  persist. *  Support  for  this  theory  comes  mainly  t rom  ballistic 


*  Finson's  theory  provides  some  targets  for  the  present  research  hut  is  neither 
a  precise  nor  a  unique  explanation  of  the  "velocity  plateau". 
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range  photographs  which  show  clearly  the  predominance  of  the  small  scale  in  the 
post-BLT  condition. 

Some  features  of  the  theory  which  the  present  research  is  in  a  position  to  prove 
or  disprove  are  as  follows: 

(a)  Velocity  "plateaus"  occur  only  paat  boundary- layer  transition. 

(b)  The  dominant  turbulence  in  the  near  wake  comes  from  the  vehicle 
boundary- layer  (provided  that  the  Re^  is  not  greatly  decreased  at 
the  shoulder  expansion  so  that  "anti- transition"  occurs). 

(c)  The  near-wake  turbulence  consists  of  scales  smaller  than  the  wake 
itself  would  generate. 

(d)  The  bulk  viscosity  of  the  small-scale  near-wake  is  decreased  below 
Che  turbulent  value,  although  it  is  still  much  higher  than  the 
molecular  viscosity. 

(e)  "Regular"  wake  behavior  sets  in  only  after  dissipation  or  diffusion 
of  the  small-scale  structure  (presumably  at  large  x/D). 

(f)  The  above  effects  are  accentuated  as  the  flight  Reynolds  number 
Rem  _  increases. 

(g)  The  above  effects  can  occur  at  virtually  any  Mach  number. 

6.3.3  THE  MEAN- FLOW  PROPERTIES 

a.  Wake  Draa  Variation.  As  explained  in  Section  6.2,  the  non-dimensional 
distance  x  behind  the  model  requires  knowledge  of  the  viscous  drag.  In  Appendix 
C  a  detailed  theoretical  estimate  of  the  cone  drag  is  made,  using  textbook  rules 
for  laminar  and  turbulent  skin  friction.  The  result  Is  a  relation  between  the 
"effective  drag  diameter"  (C0A)%  and  the  actual  cone  base  diameter  D.  If  the 
total  (viscous  and  inviscid)  drag  is  considered,  then 

(CdA)^/D  -  0.34  (6.14) 


However,  the  theoretical  rules  of  wake  behavior  demand  the  usage  of  the  "viscous" 
(or  "wake")  drag  diameter  which  is  estimated  in  the  Appendix  to  give 
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0.15D 


(6.15) 


in  the  turbulent  (post-BLT)  case  of  P^  ■  730  rim  Hg  and 
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0.14D 


(6.16) 


(CnA)  ^  - 

D  'w 

for  the  laminar  case. 

These  estimates  can  be  now  compared  directly  with  the  experimental  data  on  Figure 
40.  The  points  on  this  figure  are  computed  from  the  data  according  to  Eqns.  (613) 
and  (614).  Indeed,  despite  the  scatter,  the  points  clearly  tend  toward  the  vis- 
cour  drag  estimates.  As  for  the  initially  higher  values  of  (CqA)^,  it  is  likely 
that  this  is  simply  an  Initial  misjudgement  of  the  integration  limits  in  Eqns. 

(613)  and  (614). 

b.  Wake  Centerline  Variation.  The  wake  centerline  velocity  is  given  by  w 
defined  in  Eqn.  (612),  and  the  data  on  this  quantity  are  shown  on  Figure  41.  At 
the  lowest  pressure  (F  experiment,  P0  ■  280  mm  Hg)  for  which  the  boundary- layer 

on  the  model  is  laminar  but  where  transition  occurs  at  x/D  ~  3,  w(x)  has  definitely 
began  following  the  far-wake  data  at  x  a  30  (i.e.,  at  6  cone  base  diameter,  which 
is  about  3D  downstream  of  wake  transition).  As  one  proceeds  towards  the  post-BLT 
condition  (towards  the  A  conditions)  w(x)  departs  from  the  far-wake  behavior  sys¬ 
tematically.  In  the  A  condition  the  w(x)  at  x  ■  35  is  nearly  twice  as  large  as 
in  the  pre-BLT  condition.  Furthermore,  if  one  notes  carefully  the  WEM-A,  WEM-B 
and  WEM-C  data  (top  three  plots  on  Figure  41)  one  can  detect  a  continuing  increase, 
of  w  as  the  pressure  increases.  Therefore  it  appears  that  an  increase  in  Rei)(OJ 
increases  w  (at  fixed  x)  first  because  it  causes  transition  in  the  body  layer,  and 
subsequently  by  an  additional  mechamisin.  One  may  say  that  the  "more  mature"  the 
turbulent  boundary  layer,  the  higher  the  w(x) .  Also,  note  that  w  is  unity  at  "x  - 
3  and  decreases  abruptly  shortly  afterwards.  This  happens  because  of  the  drag¬ 
swallowing  immediately  after  the  neck.  The  above  conclusions  verify  the  theoreti¬ 
cal  points  (a)  and  (f)  (Subsection  6.3.2).  Of  course,  a  persistent,  flat  "plateau" 
is  not  observed,  perhaps  because  the  Reynolds  number  is  not  high  enough. 

c.  Viscosity  Variation.  The  theory  predicts  that  in  the  post-BLT  condition 
the  fluid  viscosity  is  reduced  considerably  below  the  turbulent  level,  although 
it  remains  equally  superior  to  the  laminar  (molecular)  level.  Figure  41  already 
constitutes  one  indirect  piece  of  evidence  that  this  is  so,  by  showing  the  expected 
increase  in  w(x)  in  the  post-BLT  case.  Furthermore,  according  to  the  far-wake 
similarity  analysis  (e.g,,  Reference  23)  the  Velocity  defect  w(x)  can  he  expressed 
in  Lerma  of  the  eddy  viscosity  vT  follows: 

w(x)  -  — (K^/ 6) 2 ^ 3  x  (6.17) 

TT 


where  R.^,  turbulent  Reynolds  number  (u^  -  u(0)  )L/v,j,. 

Since  (ue  -  u(0))  and  I.  can  be  expressed  in  terms  of  ((’|,A)’  (.n-v  Eqns.  (Oil)  and 
(6.7)),  it  follows  from  Eqn.  (6.17)  that 
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FIGURE  40.  MEASURED  DRAG  OF  TEE  HEAR 


FIGURE  41.  VARIATION  OF  HAKE  VELOCITY  DEFECT  ON  NEAR-HAKE  AXIS  WITH 
Re_  SOLID  CURVE  IS  FRGM  FAR-WAKE  EXPERIMENT. 


vT  -  const.  ~  (C_A)  ^ 

i.  D  VI 

wx 


(6.18) 


A  calculation  of  has  been  made  using  the  experimental  data  of  w,  x  and  (Cj)A)w  , 

the  latter  as  shown  on  Figure  40,  The  results  are  shown  on  Figure  42.  It  is 
abundantly  clear  that  the  bulk  kinematic  viscosity  “T  decreases  substantially  from 
the  pre-BLT  to  the  post-BI.T  case,  in  the  farther  end  of  the  wake.  Since  the  pre- 
BLT  wake  becomes  turbulent  (and  presumably  coarse-scaled)  beyond  about  3D,  this 
Implies  that  the  presence  of  shed  boundary- layer  turbulence  has  produced  the  de¬ 
creased  viscosity,  as  theory  would  have  it. 

Figure  42  shows  that  whereas  the  "regular"  vj  is  of  the  order  150-200  c.g.s.  units, 
the  decreased  level  due  to  the  small-scale  structure  is  of  order  50-100  c.g.s. 
units.  A  rough  calculation  of  the  expected  laminar  value  gives  1-10  c.g.s.  units. 
Thus  the  small-scale  turbulence  brings  down  but  still  keeps  it  substantially 
higher  than  the  laminar  level,  finally,  the  relation  among  these  various  levels 
can  be  checked  by  additional  existing  data,  From  the  computation  performed  of 
the  "wake  Reynolds  number"  Re  •  (u#  -  u(0))L/ve  (see  Eqn.  (6.10))  we  typically 
find  Rew  ~  1000  for  the  pre-BLT  case.  ThuB  the  ratio  o  molecular- to-eddy  vis¬ 
cosities  is 


Re 


w 


15 

-  750  - 


_1 

50 


(6.19) 


which  compares  well  with  the  number  quoted  above  in  c.g.s.  units. 

d.  The  Turbulent  Prandtl  Number.  Although  not  of  direct  importance  to  the 
theory  insofar  as  the  type  of  wake  structure  is  involved,  the  Prandtl  number  crT 
controls  the  predicted  fluid  temperature  8  in  the  wake  through  the  familiar  formula 


ctt(y  *  1)  Me  w 


(6.20) 


The  crip  has  been  computed  from  the  data  and  is  shown  on  Figure  43.  It  tends 
asymptotically  in  all  cases  to  the  well-established  level  of  about  0.85  charac¬ 
teristic  of  axisymmetric  flow  (Reference  23). 

e.  Radial  Profiles  of  Velocity  and  Temperature.  Eventually,  any  predictive 
theory  of  the  near  wake  will  either  assume  or  predict  the  radial  velocity  u1  and 
temperature  T  variations.  It  is  usually  expected  that  the  post-BLT  wake  will  have 
rather  novel  variations  of  thia  type.  The  data,  as  computed  from  Eqn.  (b.6)  and 
presented  on  Figures  44  and  45,  however,  show  that  u  and  X  are,  in  fact,  quite 
"Gaussian".  The  solid  lines  shown  are  the  self-similar  variations 
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-0.34  Tl2 
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which  were  found  by  Philco-Ford  in  the  far  wake  (Reference  23).  It  is  quite  sur¬ 
prising  that  the  data  actually  reaches  this  behavior  almost  immediately  after  the 
"neck"  (e.g.,  at  x  ■  2D);  the  "bulge"  at  T)  >  2  is  due  to  the  shed  boundary-layer. 

6.4  TURBULENCE  PROPERTIES 

6.4.1  METHOD  OF  COLLECTING  THE  DATA 

a,  Principle.  We  are  interested  in  mapping,  first,  the  axial  and  radial 
distribution  of  the  velocity  and  temperature/density  fluctuations,  and  secondly, 
the  distributions  of  the  integral  scales  of  the  longitudinal  velocity  and  temper¬ 
ature  fluctuations.  The  principle  of  measuring  the  fluctuation  intensities  is  as 
follows.  The  hot-wire  is  positioned  at  the  desired  point  in  the  wake  and  its  wide¬ 
band  rms  output  e  and  mean  voltage  drop  V  are  measured  at  each  of  a  number  of 
electrical  heating  currents  i.  The  shape  factors  of  the  curve  defined  by  the  plot 
e/V  versus  i  gives  the  desired  fluctuation  magnitudes.  For  the  scale  measurement, 
one  invokes  Taylor's  hypothesis  that  the  scale  can  be  extracted  from  the  measure¬ 
ment  of  the  power  spectrum  e2(f)i 

A  ~  (6.22) 

e 


2 

where  e  (f  -  0)  is  the  spectral  density  at  zero  (or  very  low)  frequency  and  e  is 
the  wideband  fluctuation  already  measured  above.  Therefore  one  can  simply  measure 
e2(f  ■  0)  at  that  point  and  combine  it  with  the  e  already  measured  to  come  up  with 
the  scale  length.  We  shalL  now  give  a  somewhat  more  detailed  description  of  how 
this  is  done  in  the  laboratory. 


b.  Description  of  Measurement  (Refer  to  Figure  46).  At  each  pressure  a 
number  of  positions  (x  stations)  are  chosen  for  study.  Normally  these  positions 

are  in  increments  of  0.2D  starting  at  x  -  0  (that  is  0,  0.2,  0,4, . 1.4,  1.6) 

up  to  x  *  1.6D;  then,  in  increments  of  0.4D  from  2D  to  6D  (i.e,,  2,  2.4, . 5.6, 

6).  At  each  x  data  are  then  taken  in  the  radial  direction  at  points  spaced  0.005D 
apart,  beginning  with  point  yQ  and  ending  with  point  y  ,  Typically  for  x/D  >  1.4, 
y  -  0  (i.e.,  the  first  point  is  read  on  the  wake  axis)  and  ye  ■  0,395  D.  Nearer 
tne  body  y0  and  ye  are  chosen  to  bracket  the  active  portions  of  the  flow.  For  ex¬ 
ample,  at  x  *  0  (i.e,,  at  the  plane  of  the  cone  base)  the  measurement  starts  at  yo 
■  0.525D  (i.e.,  0.025D  above  the  cone  surface)  and  continues  to  ye  *  0.695D,  i.e., 
0.195D  above  the  cone  surface. 

Data  are  punched  on  IBM  cards  in  real  time.  The  hot-wire  probe  current  i  is  first 
fixed  at  the  desired  value.  As  the  wind-tunnel  runs  the  probe  is  first  located  at 
the  desired  starting  position  (say  at  x  -  2D,  y/D  *  0)  and  then  a  switch  is  thrown 
which  begins  traversing  the  probe  in  the  radial  position  at  constant  speed,  Every 
0.0Q5Dof  probe  travel  a  pulse  automatically  triggers  the  data-acquisition  system 
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and  a  pair  of  numbers  1b  recorded  on  IBM  cards.  At  the  end  of  travel  (that  is, 
when  the  data  at  the  end  point  ye  have  been  punched)  the  probe  is  returned  to  its 
starting  point  y0,  the  wire  heating  current  i  changed  to  the  next  desired  level, 
and  the  traverse  is  repeated.  When  the  same  path  has  been  traversed  for  all  the 
desired  currents,  the  probe  is  moved  to  another  axial  x-station,  and  the  process 
is  repeated. 

The  Philco-Ford  A/D  conversion  system  presently  has  two  active  channels,  meaning 
that  at  each  point  (x,y)  only  two  numbers  can  be  recorded.  However,  to  get  both 
the  fluctuation  intensities  and  the  scales  we  saw  in  the  preceding  section  that 
we  need  three  measurements,  i.e.,  e,  V  and  e2(f  =0).  To  resolve  this  problem, 
it  was  decided  to  multiplex  the  data  acquisition  so  that  at  one  point  one  reads 
e  and  V,  whereas  at  the  next,  one  reads  e,  e2(f  *  0),  then  back  again  to  e,  V, 
etc.  A  special  multiplexer  was  built  to  perform  this  switching  function. 

Referring,  now,  to  Figure  46,  the  hot-wire  signal  is  channeled  to  an  ADP  12/13 
Hot-Wire  Amplifier,  which  boosts  it  to  a  gain  G  of  100-500,  compensates  it  for 
thermal  lag  and  filters  its  components  below  3  Hz  and  above  1  MHz.  The  current 
1  of  the  wire  is  permanently  displayed  during  the  run  on  a  digital  voltmeter; 
the  mean  wire  voltage  V  is  channeled  to  the  multiplexer,  as  well  as  to  a  flow- 
calibration  analog  computer  (Model  ADP-17).  During  the  earlier  part  of  the  work 
the  latter  was  used  to  produce  flow  calibrations  of  the  anemometer  prior  to  each 
measurement  (see  further  discussion  of  this  under  Section  III. 

The  a.c.  hot-wire  signal  e  is  carried  to  a  Hewlett-Packard  Model  3400A  rms  volt¬ 
meter,  through  which  it  is  converted  to  a  d.c.  voltage  Z  with  a  transfer  func¬ 
tion  tj_  such  that  e(rms  volts)  **  t^Z(d.c.  volts).  Usually  the  rms  meter  is  set 
to  a  3-V  full-scale  so  that  t.  *  3.  The  output  Z  is  then  fed  to  a  Hewlett-Packard 
Model  3440A  Digital  Voltmeter  ( DVM)  which  is  part  of  the  Philco-Ford  A/D  system. 

In  the  meantime  the  a.c.  wire  signal  e  is  also  channeled  to  a  Hewlett-Packard 
model  302A  Wave  Analyzer  tuned  to  a  frequency  f  which  usually  lies  between  100 
and  250  Hz.  It  is  known  that  the  signals  involved  are  very  high  in  frequency, 
so  that  100  <  f0  <  250  Hz  lies  well  into  the  low-frequency  "flat"  portion  of  the 
spectrum.  Thus  f0  can  be  considered  to  give  the  zero-frequency  end  e2(f  ■  0)  of 
the  spectrum.  The  function  of  the  wave  analyzer  ia  to  filter  the  wideband  signal 
e  with  a  passband  of  7  Hz  around  fQ  (i.e.,  3.5  Hz  on  either  side  of  f  ,  with  a 
-3  db  drop).  The  signal  passing  through  this  window  is  then  extracted  from  the 
analyzer  in  the  form  of  an  a.c.  signal  eQ  with  a  transfer  function  J  such  that 


eQ(folts  rms)  -  J  ei(volts  rms)  (6.23) 

The  transfer  function  J  is  usually  set  at  33.3,  so  that  a  30  mV  equivalent  (sine 
wave)  input  e^  will  produce  a  one-volt  eQ.  If  e2(f)  is  the  spectral  density 
(i.e.,  the  integral  of  e2(f)  over  f  produces  the  wideband  signal  e~)  then  the 
equivalent  input  e^2  is  7e2(f  ).  Therefore  the  desired  quantity  e2(fQ)  is  re¬ 
lated  to  the  analyzer  output  e02  by 

2  e  2 

e  if  )  "  (6.24) 
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Next,  the  analyzer  output  e0(rtns  volts)  Is  fed  to  another  Hewlett-Packard  Model 
3400A  rms  meter  with  transfer  function  t?  (usually  t2  ■  3)  such  that  its  d.c. 
output 

t»F(volts  d.c.)  -  e  (volts  a.c.)  (6.25) 

l  o 


Thus  collecting  the  various  transfer  functions  in  the  sequence,  the  desired 
spectral  density  e^(f  )  is  found  from  the  reading  (or  F)  by 


(6.26) 


The  output  F  is  also  channeled  to  the  multiplexer.  This  concludes  the  description 
of  essential  elements  in  the  measuring  system  shown  on  Figure  46. 

As  already  mentioned,  each  traverse  begins  with  the  probe  poiBed  Just  before  the 
first  y  position  (i.e.,  y  )  and  the  multiplexer  set  to  pass  through  the  V  (wire 
mean  voltage)  signal.  Thus  as  soon  as  probe  motion  begins,  the  first  two  numbers 
punched  are  V  and  Z,  i.e.,  the  data  needed  to  get  the  fluctuation  magnitudes. 

The.  multiplexer  next  switches  to  produce  the  F  signal,  so  that  at  the  second 
point  Z,  F  are  read.  The  multiplexer  alternate  switching  action  continues  to  the 
end  of  the  traverse  so  that  one-half  of  the  points  read  are  Z,  V  and  the  other 
half  are  Z,  F  data.  Thus  the  fluctuation  magnitudes  are  obtained  (from  Z,  V)  at 
points  spaced  2  x  0.005D  -  0.010D  apart  and  so  are  the  scales  (from  Z,  F) ,  the 
two  alternating  along  the  traverse.  Typically,  80  points  are  probed  along  each 
radius,  giving  40-point  maps  for  each  of  the  magnitudes  and  the  scales. 

c.  Data  Reduction  Procedure.  The  turbulence  data  are  reduced  with  the  aid 
of  computer  program  WEM-1V  (DRB4E) .  The  following  are  notes  explanatory  to  the 
program,  with  the  sub-heading  numbers  referring  to  the.  similarly  numbered  program 
steps  (see  Appendix  D) . 

(7.3.1)  The  "current  group"  refers  to  the  wire  heating  current,  which 
of  course  changes  from  one  wire  to  another. 

(7.3.3a)  The  probes  used  in  the  experiment  have  a  line  resistance  of 
1.42  ohm. 

(7.3.4)  For  low  overheats,  it  is  known  that  R  is  a  linear  function  of 
i.2 ;  R  -  c  +  c^i  .  However,  for  higher  overheats  the  varia¬ 
tion  is  quadratic,  which  explains  the  presence  of  the  term. 

(7.3.11)  The  "Reynolds  No."  quantity  printed  in  the  output  of  the  WFM-II 
program  is  based  on  a  hot-wire  diameter  of  0.00001"  (ten  micro¬ 
inches).  The  wire  used  in  the  present  (WEM-IV)  tests  are  of 
0.00002"  diameter,  which  accounts  for  the  factor  of  2  in  this 
step. 


-88- 


(7.3.11a/  Originally,  the  hot-wire  was  "flow  calibrated"  in  the  stream 
7. 3. lie)  external  to  the  wake  so  that  its  Nusaeit  number  Nu  dependence 
on  the  local  Reynolds  number  Re  could  be  learned.  Then,  once 
this  "flow  calibration"  Nu  ■  Nu(Re)  was  known  the  Reynolds 
number  from  step  (7.3.11)  was  used  to  find  Nu  at  every  point 
of  measurement;  from  this  Nu  and  the  Nu  ~  Nu(Re)  relation  the 
quantities  (3Nu/SRe) (Re/Nu) ,  etc.,  were  computed  for  use  in  the 
sensitivity  coefficients  (Reference  17).  However  it  was  found 
that  the  radial  Nu  profile  so  computed  differed  from  that  found 
by  a  direct  measurement  Num  of  the  Nuaselt  number  in  the  wake: 


Nu 
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where  the  quantities  on  the  right  are  actually  measured  at  each 
point.  It  was  furthermore  found  that  if  Nu  was  used  instead 
of  Num,  the  turbulence  results  were  suspect  and  varied  greatly 
from  one  wire  to  another  (for  example,  the  cross-correlation 
coefficient  rQT  was  found  to  be  positive,  contrary  to  all  pre¬ 
vious  experience  and  all  expectations).  The  situation  was 
rectified  by  discarding  the  flow-calibration  process  entirely 
in  favor  of  using  the  Num  matched  to  the  local  Re  from  the  WEM-II 
program  results.  Steps  (7.3.11a)  to  (7. 3. lie)  in  the  program, 
therefore,  aim  at  deriving  the  values  of  the  Nusselt  number  and 
recovery  factor  from  the  known  hot-wire  data  at  each  point.  It 
will  be  noted  that  the  curve-fit  for  Num  *  Num(Re)  is  linear 
(step  7.3.11c),  since  such  a  linear  fit  was  Judged  best  from 
the  available  data,  (over  the  range  of  the  Reynolds  numbers 
used).  The  recovery  factor  71  s  Taw/T0  is  fitted  with  a  quad¬ 
ratic  in  Re  also  as  a  result  of  inspection  of  the  available 
data. 

6.4.2  PRELIMINARY  RESULTS  OF  THE  TURBULENCE  MEASUREMENTS 

A  portion  of  the  data  matrix  has  already  been  covered,  and  preliminary  results 
from  such  measurements  will  be  discussed  in  this  Section.  It  should  be  noted 
that  the  matrix  requires  "mass-production"  of  quantitative  hot-wire  anemometer 


*  In  most  other  experiments  conducted  at  Philco-Ford  the  Nu-Re  curve-fit  is  of 
the  type  Nu  *  A  +  B  /"Re  +  C  Re.  The  present  linear  curve-fit  (which  essen¬ 
tially  sets  B  ■  0)  reflects  only  the  present  Re  range  and  should  not  be  mis¬ 
taken  as  a  more  generaL  hea t- transfer  correlation.  Also,  the  recovery  factor 
is  presently  fitted  by  7)  ■  +  H2Re  +  H^Re^  because  it  attains  a  maximum 

near  Re  ■  1.  This  maximum  is  due  to  the  finite  wire  length,  by  which  7]  does 
not  anymore  keep  increasing  fs  Re  is  lowered  because  end  losses  take  over. 
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data  In  an  unprecedented  scale,  The  data  are  therefore  to  be  considered  tenta¬ 
tive  since  a  need  for  improving  the  technique  has  been  found  to  exist  in  many 
areas.  The  discussion  will  center  on  Figures  47  through  49. 

The  first  important  finding  concerns  the  fluctuation  magnitudes  u'/u  and  T'/T 
in  the  cone  boundary-layer  at  *  «  0,  i.e,,  at  the  soulder  (see  Figure  47).  The 
velocity  fluctuation  u'/u  decreases  monotonical iy  from  the  wall  outwards,  while 
T'/T  has  a  peak  in  the  middle  of  the  boundary- layer ,  reaching  T'/Ttr  0.09.  This 
behavior  is  consistent  with  expectations  from  Kistler's  work  (Reference  9)  if  one 
assumes  no  Reynolds  number  or  geometry  effect,  That  is,  the  present  data  at  Rq 
■  1.170  and  a  cone  geometry  agree  numerically  with  Kistler's  data  for  a  flat 
plate  taken  with  a  much  larger  Rq  (  *  20,000).  Thus  the  near  wake  of  this  test 
begins  "normally"  at  the  body  shoulder.* 

The  second  noteworthy  result  is  that  the  velocity  fluctuations  u'/u  decay  abruptly 
beyond  the  shoulder,  consistently  with  the  theoretical  prediction  on  the  effect 
of  the  expansion  on  the  shed  turbulence;  that  is,  from  a  maximum  level  of  about 
5  to  67.  (u'/U  a  0.06)  in  the  boundary  layer,  these  fluctuations  decrease  to  the 
1-27.  level.  At  least  to  a  certain  degree  this  represents  the  "anti-transition" 
(re-laminarization)  effect  by  which  turbulence  vanishes  upon  a  Budden  Reynolds 
number  decrease.  Now,  Finson  (Reference  18)  remarks  that  a  minimum  turbulence 
Reynolds  number  Re/\  of  about  35  is  needed  to  guarantee  survival  of  turbulence  in 
the  near  wake;  this  threshold  number  agrees  with  a  general,  rough  estimate  made 
in  Appendix  E  of  this  report.  The  numerical  data  quoted  above  imply  that  some 
turbulence  activity  has  survived  in  our  near  wake.  Thus  the  Re/\  in  our  case 
(Figure  47)  must  be  at  least  35  or  so  in  the  near  wake,  and  presumably  much  higher 
in  the  boundary- layer  (at  x  ■  0) .  Indeed,  using  u'  from  Figure  47  (the  x  *  0 
case),  the  corresponding  mean-velocity  data  (from  Reference  19)  and  an  estimated 
A  •  6/5  we  compute  Re^  ■  150  in  the  boundary- layer .  The  conclusion,  again,  is 
that  for  the  A  experiment  of  our  tests,  the  configuration  fulfills  the  "small- 
scale  structure"  requirements  envisaged  by  the  theory. 

In  Figure  47  it  is  also  seen  that  beyond  the  shoulder  (x  >  0)  the  velocity  fluc¬ 
tuations  become  very  large  toward  the  bottom  of  each  graph,  as  soon  as  the  "re¬ 
circulation  region"  is  entered.  Two  comments  are  in  order  here.  First,  in  the 
latter  region  the  measurements  cannot  be  accurate  since  the  direction  of  fluid 
motion  is  quite  irregular.  Secondly,  the  fluctuations  are  plotted  in  terms  of 
u'/u  where  "u  is  the  local  mean  velocity;  the  latter  is  very  small  in  that  region. 
That  is,  u'  is  not  necessarily  large  in  itself.  The  computer  program  presenting 
the  data  of  Figure  47  also  prints  the  fluctuations  u'  and  T'  normalized  with 
other  quantities  such  as  the  flight  speed  ue,  etc,  The  normalizations  u'/u^, 
etc.,  will  be  discussed  in  a  later  report. 


*  The  point  made  here  is  that  the  wake  initial  conditions  begin  with  a  "familiar" 
profile.  Since  both  this  test  and  Kistler's  data  are  at  low  Me  (  ~  3),  theBe 
results  leave  unanswered  the  Mach  number  dependence  of  the  wake  initial  condi¬ 
tions.  Such  an  answer  is  forthcoming  from  the  JPL  and  AEDC  tests  (see  Sections 
II  and  III).  Also,  in  Section  IV  we  summarize  the  stnte-of-knowledge  based  on 
such  results. 


90- 


°o0oa<t>e'JO°°°°' 


O  u'  •  T 


WIDEBAND  TEMPERATURE  AND  LONGITUDINAL  VELOCITY  FLUCTUATIONS 
IN  THE  NEAR- WAKE. 


FIGURE  ^9.  VELOCITY  FLUCTUATIONS  IN  POST-BLT  NEAR  WAKE  AS  VIEWED  IN 
FLUID  AND  GROUND-FIXED  COORDINATES. 


The  temperature  fluctuation  behavior  shown  on  Figure  47  1b  substantially  differ¬ 
ent  from  that  of  u'/u.  The  fluctuations  T'/T  are  unaffected  by  the  expansion  and 
continue  towards  the  wake  neck,  at  x/D  ~  1,  with  intensity  undimi nished  at  nearly 
10%  maximum.  This  result  confirms  the  now-familiar  ballistic-range  observations 
showing  the  fine-grained  boundary- layer  structure  "pouring"  into  the  wake,* 

Thus,  subject  to  the  Re^  limitations  of  this  test,  it  appears  that  a  substantial 
amount  of  boundary- layer  turbulence  continues  on  into  the  wake,  but  that  it  con¬ 
sists  chiefly  of  temperaturu/density  fluctuations.  Qualitatively,  this  result 
could  have  been  predicted  in  advance  because  of  the  well-estabLished  "Strong  Rey¬ 
nolds  Analogy"  connecting  u'  and  T'* 


Ll 

T 


(V  -  1)M2  ^ 


(6.28) 


Since  the_local  Mach  number  M  increases  upon  expansion,  the  temperature  fluctua¬ 
tions  T'/T  increase  in  a  way  tending  to  balance  the  decreasing  u'/u. 

In  addition  to  the  peak  in  T'/T  of  Figure  47  which  was  just  above  identified  with 
the  boundary- layer  remnant,  a  aecond  peak  appearB  between  the  remnant  and  the  wake 
axis.  This  peak  is  intense  and  narrow  and  is  obviously  associated  with  the  free- 
shear  layer.  Its  prominence  could  perhaps  be  expected,  but  the  absence  of  equally 
intense  velocity  fluctuations  at  the  same  location  is  puzzling,  especially  since 
u'/u  does  not  appear  like  a  well-defined  "peak"  there. 

The  evolution  of  the  fluctuations  u'/u  and  T'/T  beyond  the  neck  is  shown  on  Fig¬ 
ure  48.  Just  beyond  the  neck  (at  x/D  "  1.4)  the  boundary- layer  remnant  and  the 
shear  layer  can  be  distinguished  from  each  other.  Subsequently  they  merge  into 
each  other,  the  merging  process  being  complete  about  2D  beyond  the  neck  (at  x  - 
3. 2D).  Note  that  there  appears,  however,  to  remain  a  "velocity  fluctuation  wake" 
distinguishable  from  a  "temperature  fluctuation  wake"  at  x  ■  3. 2D.  Outwardly, 
this  is  evidence  of  a  "double"  structure  wake  whereby  strong  velocity  fluctuations 
at  the  v_cre  are  surrounded  by  an  annulus  of  temperature  fluctuations^  In  turn, 
thie  may  be  taken  as  evidence  that  the  wake  is  formed  by  (a)  the  T'/T  from  the 
remnant  (which,  as  we  saw  earlier,  survives  the  shoulder  expansion)  and  (b)  the  u'/u 
contribution  from  the  shear  layer,  However,  judgement  on  this  conclusion  should 
be  withheld  since  the  shear  layer  shows  no  clear  prominence  of  u'/u  and  because 
the  shear  layer  already  has  large  T'/T  (observe  the  first  graph  of  Figure  48) 
before  it  mixes  with  the  remnant.  In  fact,  the  double  structure  of  the  wake  at 
x  *  3. 2D  is  a  simple  consequence  of  adiabatic  mixing  of  compreosible  turbulence 
as  can  be  seen  by  re-writing  Eqn.  (6.28)  as  follows: 


*  The  fact  is  often  overlooked  that  the  ballistic-range  photographs  show  the 
development  of  density  and  temperature  fluctuations,  not  of  velocity  fluctua¬ 
tions,  in  the  near  wake.  The  interpretation  of  (velocity)  turbulence  behavior 
solely  by  means  of  such  photographs  is  therefore  quite  out  of  order. 


t 
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-  (V  -  1)M2  -  0.4  M2  <  1  for  M  <  1.56 

u'/u 


Since  the  axis  streamline  begins  at  the  neck  with  M  -  0  and  Increases  thereafter, 
the  velocity  fluctuations  on  the  axis  will  be  greater  than  the  temperature  fluc¬ 
tuations  for  some  distance  past  the  neck,  i.e,,  until  M  (axis)  ■  1.5,8.  Further¬ 
more,  since  M  is  higher  outside  the  wake  the  two  curves  u'/u  and  T ' / T  should 
cross,  which  is  shown  by  Figure  48  to  occur  in  practice, 

Figure  49  eummariees  the  tentative  fluctuation  intensity  u'/TT  results.  Note  the 
relative  roles  and  intensities  of  the  boundary- layer  remnant  and  of  the  shear 
layer.  Note,  also,  that  a  fluctuation  lntenaity  peak  occurs  in  the  neighborhood 
of  x/D  «  2 . 5D,  reaching  97.  of  the  local  velocity  u,  and  6%  of  the  external  velo¬ 
city  ue.  This  peak  has  also  been  observed  by  Oran  (Reference  24).  However, 
Gran's  observation  that  the  fluctuation  maximum  occurs  on  the  axis  is  not  borne 
out  by  the  present  results. 
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SECTION  VII 


ELECTRON  SPACE-TIME  CORRELATIONS  IN  TURBULENT  PLASMAS 


In  the  present  reporting  period  the  study  of  the  relation  between  gas  temperature 
and  eLectron  density  fluctuations  in  a  free  turbulent  flow  has  been  completed. 

Most  of  the  effort  concentrated  on  assembling  the  data  collected  in  previous  con¬ 
tract  periods,  reducing  supplementary  data  when  necessary,  and  collecting  the  re¬ 
sults  into  a  technical  report.  This  report  is  now  under  preparation  and  will  be 
issued  soon. 

The  major  findings  of  this  study,  abstracted  from  the  report  just  mentioned,  are; 

(1)  The  correlated  eddies  of  the  gas  temperature  tend  to  be  near- spherical 
structures  throughout  the  turbulent  volume,  with  a  scale  length  of 
the  order  of  1/5  the  flow  diameter. 

(2)  By  contrast,  the  electron  eddies  are  thin  structures  elongated  in  the 
flow  direction,  especially  on  the  axis.  The  "aspect  ratio"  of  these 
elongated  eddies  is  on  the  order  of  1:5.  Near  the  flow  edge,  however, 
the  elongation  diminishes  and  the  electron  eddies  appear  more  coinci¬ 
dent  with  and  of  similar  shape  as  the  temperature  eddies. 

(.3)  The  electrons  decay  to  a  considerable  degree  within  a  distance  of  about 
1/2  the  flow  diameter,  that  is,  recognizeable  lifetime  of  an  electron 
eddy  is  of  order  d/2. 

The  above  results  are  specific  to  the  particular  experiment  performed,  which  is 
described  in  detail  in  References  25,  26  and  27.  A  very  detailed  picture  of 
simultaneous  gas- dynamical,  thermal  and  electronic  behavior  thus  exists  for  a 
high-temperature  turbulent  flow.  For  example,  simultaneous  temperature  and 
electron  density  signals,  recorded  at  many  axial  and  radial  positions  in  the  flow, 
exist  on  a  number  of  1/2-inch,  10, 000-foot  tapes.  This  tape  library  is  avaiLabLe 
for  further  data  processing  without  the  need  to  set-up  and  repeat  the  tests. 

The  actual  chemistry  of  the  flowfield  is  perhaps  the  most  difficult  re-entry 
feature  to  duplicate  in  a  ground  facility,  and  the  results  obtained  can  only  be 
applied  to  the  flight  case  via  a  theoretical  calculation.  That  is,  if  such  a 
theory  is  available  for  predicting  electron  fluctuations  in  flight,  it  can  use 
the  results  of  this  section  as  a  validity  test.  However  no  such  theory  is  pre¬ 
sently  available.  During  the  contract  period  just  ended,  an  effort  was  made  to 
prepare  such  a  theory  for  non-equilibrium  reacting  mixtures  along  lines  similar 
to  those  of  Reference  28  for  equilibrium  mixtures.  Progress  was  indeed  made  as 
far  as  to  indicate  the  faults  of  previous  attempts  (References  29  and  30  for 
example)  and  some  major  improvements.  Calculations  actually  began;  however  the 
effort  to  complete  them  was,  as  expected,  much  beyond  the  resources,  and  out¬ 
side  the  work  statement,  of  this  program.  The  theoretical  work  was  therefore 
temporarily  halted. 
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APPENDIX  A 


EFFECT  OF  LATERAL  DISPLACEMENT  OF  PROBES  ON  TRUE  PROBE 
HEIGHT  ABOVE  CONE  SURFACE 


The  four  probes  used  in  the  PF/AEDC  cone  boundary- layer  measurements  consisted 
of  a  hot-wire  anemometer  (W),  a  pitot  probe  (P),  a  thermocouple  (T)  and  a  sur¬ 
face  sensing  probe  (S).  These  probes  were  mounted  on  a  single  rake  with  their 
sensing  tips  located  in  a  plane  normal  to  the  cone  surface  (we  can  ignore  small 
displacements  in  the  atreamwiBe  direction).  In  order  to  avoid  contact  between 
the  measuring  probes  V),  F  and  T  and  the  cone  surface,  these  probes  were  displaced 
vertically  above  the  tense  probe  (l.e.,  in  the  direction  normal  to  the  cone  sur¬ 
face)  by  an  amount  y0(S),  y  (P) ,  and  y  (T)  respectively  (the  distance  yQ(i)  la 
the  separation  between  the  bottom  of  tne  surface  probe  and  the  center  of  the 
sensing  area  of  the  measuring  probe).  Thia  is  shown  schematically  in  Figure  A1 
where  the  probes  have  been  located  with  the  surface  sensor  contacting  the  cone, 
Figure  Al.  also  indicates  the  lateral  displacement  of  the  probes, 

Because  of  surface  curvature  the  true  height  H  of  any  probe  above  the  cone  sur¬ 
face  will  differ  from  the  aun  y  +  y0 (where  y  is  the  vertical  travel  of  the  probe) 
and  is  given  instead  by  the  quantity  H  ■  (R  -  Rc)  where  R  is  the  length  of  the 
straight  line  segment  between  the  cone  center  and  the  probe  and  RQ  is  the  cone 
radius.  This  is  illustrated  in  Figure  A2  where  the  probe  position  has  been  de¬ 
fined  by  the  polar  coordinates  R  and  6, 


The  true  probe  height  can  be  found  in  the  following  manner.  From  Figures  A1  and 
A2  it  can  be  shown  that: 


Since: 


then: 


and; 
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so  that: 


H(W) 

■  R(W)  -  R 

c 

H(P) 

-  R(P)  -  R 

c 

(A- 3) 

H(T) 

-  R(T)  -  R 

It  is  clear  that  the  effect  of  surface  curvature  (i.e.,  the  difference  between 
(y  +  y0)  and  H  increases  with  decreasing  K-Station  and  with  decreasing  y. 


Since  the  angles  0(i)  ate  very  amall  H(i)  can  be  approximated  by  the  quantity 
D(i)  which  ia  defined  in  Figure  A3  where,  for  simplicity,  only  the  pitot  probe  is 
represented.  It  can  be  shown  that! 

8(P)_  -  L(P)/R 
o  c 


flnd 

so  that: 


R0Cp) 


L(P)/ tan  0(P) 


D(P) 

-  R0  ■  Vp)  +  vP> 

+  y 

that : 

D(W) 

-  Rq  -  Ro(W)  +  yQ(W) 

+  y 

D(T) 

-  Rq  -  Rq(T)  +  yo(T) 

+  y 

(A-4) 


These  expressions  have  the  advantage  that  R  -  R  (P>,  the  correction  for  surface 
curvature,  is  independent  of  y.  o  o 


An  estimate  of  the  effect  of  surface  curvature  ia  provided  by  considering  the  x  - 
19  inch  station. 


Let : 

R 

c 

1.75  inches 

y0(w)  - 

0.030  inches 

L(S)  - 

i(S)  -  0.0625 

inches  yQ(P)  ■ 

0.010  inches 

L(P)  - 

L(T)  -  0.1875 

inches  yo(T)  ■ 

0.040  inches 

corresponding  values  of  y(i),  H(i)  and  D(i)  are  listed  below: 

y(i)  - 

0 
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y(i) 

H£ii 

D(l) 

Hji hr  i) 

P ill  -  H( i) 

W 

.030 

.0299813 
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-1.87  x  lO”5 
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P 
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.0188493 

.0189572 

.0088 

.000108 

T 

.040 
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.0489572 

.00868 

.000274 
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FIGURE  A3.  APPROXIMATION  FOR  TRUE  PROBE  HEIGHT  ABOVE  CONE  SURFACE. 


y(i)  - 

0.25 
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It  is  apparent  that  the  surface  curvature  effect  is  most  significant  for  the  pitot 
probe  located  near  the  surface  where  the  error  H(P)  -  y(P)  Is  of  the  same  order  as 
y(P).  It  is  necessary  therefore  to  Include  the  correction  for  surface  curvature 
in  the  data  reduction  procedure.  It  is  also  clear  that  the  difference  D(l)  - 
H(i)  is  within  the  precision  with  which  the  probe  position  can  be  specified. 
Therefore  we  can  set  H  ■  D  and  use  the  expressions  (A-4)  to  evaluate  the  "true" 
probe  height. 
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APPENDIX  E 


X-PROBE  RESPONSE  WITH  PRESSURE  FLUCTUATIONS 


The  theory  of  the  X-probe  response  at  high  speeds  and  in  the  absence  of  pressure 
fluctuations  has  been  described  in  Appendix  B  of  Reference  5.  When  pressure 
fluctuations  cannot  be  neglected  we  must  recognize  that  the  mode  and  flow  variables 
are  no  longer  identical  (see  Reference  1)  but  rather  related  by  the  expressions: 


T1 

-  o' 

+  (Y  -  1)  tt' 

(B-la) 

u* 

m  T  1 

+  (nx/M)  tt' 

(B-  lb) 

v’ 

“  S' 

+  (ny/M)  tt' 

(B-lc) 

P' 

-  Yrr 

1 

(B-ld) 

For  a  wire  normal  to  flow  (i.e,,  not  sensitive  to  v')  the  response  of  the  wire  is 
given  by: 


e'  -  e,rTT'  -  emm'  (B-2) 

Expanding  T^'  and  m'  in  terms  of  T',  u',  and  p',  and  introducing  Eqns  B-la  through 
B-ld  and  the  definitions 

y  =  - 1 - -  .  (3  s  CY  “  lJM2^ 

1  +  C(Y  -  l)/2]*r 


equation  B-2  can  be  re-written  as 


e'  -  e  a'  +  e  T'  +  [><Y  -  1)  +  a-(Y  -  l)Mn  )e„ 
O'  xi 


(B-3) 


-  <*r  +  V'J  n' 


Defining  the  bracketed  term  as  the  sensitivity  to  pressure  fluctuations,  e  , 
Eqn,  (B-3)  becomes: 

e 1  ■  e  o '  +  e  T '  +  e  n ' 

a  T  tt 


where  e  is  defined  in  the  Horkovin  AGARDO-graph  (Reference  17)  and  it  can  be 
shown  that: 


e 

TT 

e 

a 


.  (Y  -  1H1T 


1  +  (Y  -  1)M 


ji  „  r'1  -  ",M) 


(M 


1  +  (Y  -  1)M 


X 
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and  X  la  Che  conventional  mode  parameter  e^e^. 


For  a  skew  wire  u*  la  replaced  with  u*  +  v'  cot  0  (see  Reference  5)  and  the  sen- 
ltl.itj  coefficient#  are  based  now  on  the  velocity  component  normal  to  the  wire , 
r^,  nod  not  the  total  velocity  vector,  u.  In  this  o...  ..  can  .rite  th.t  .gain. 


e1  -  e_T_,'  -  e  m1 

'T  T'  m 


which  after  some  algebraic  manipulation  reduces  to: 

e'  -  «  ff»  +  eTT'  +  eT  cot  0  5 '  +  [C«(V  “  D  +  <%<Y  “  DM 

+  «(v  -  l)nyM  cot  0]eT  ^  cot  0  +  1]  ej  tt ’ 

If  we  let  nx  -  -  ^  then  Eqn.  (B-4)  becomes  finally: 

e'  .  e^a'  +  eTT'  +  cot  0  ?  '  4-  eMi' 


n 

+  [ff(Y  -  l)nyM  eT  -  j*  •„]  cot  0  n* 


Let  us  define  now 


Since : 


A  a  »(V  i)«yM  eT  -  ^ 


e  ■  +  e  ( 

a  T  m 


e  ■  ® 

T  T  m 


(B-4) 


(B-6) 


— (e  +  e_)  and 

<k  +  0  a  T 


S+’V  e<3  “  B+T(  t 


Consequently,  L  can  be  written  as. 


-  of<Y  -  i>v  [r+*  ea '  eJ 


From  the  definition  of  «  and  fl  it  can  be  shown  that: 

fv  -  1)M2  +  1 

-v  +  8  -  — u - - ~2 

1  +  [(v  -  1 ) / 2 3 M 


(B-7) 


7+9  1  +  <Y  -  DM2 


*L.S-  C 
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-*  +  9 


— (V-  -  DM^ 

1  +  (Y  -  1)M2 


and  therefore  that: 


Vn  M  -  (M2  -  l)r* 

— *  . — —  e  “  h(M,n  )  e 

i  +  <y  -  dm2  T  y  T 


Substituting  Eqn.  (B-8)  into  Eqn.  (B-5)  we  have: 


e  c? '  +  e_T'  +  e  cot  0  t1  +  e  tt '  +  6  cot  0  it* 
a  T  T  5  n 


(B-8) 


(B-9) 


From  Morkovin  (Reference  17)  the  hot-wire  sensitivity  to  the  sound  mode  for  a 
wire  normal  to  the  flow  is  given  by: 


where 


e  -  f (M)  e  +  g(M,n  )  e 

TT  0  X  T 


(Y  -  1HM  -  1) 


1  +  (y  -  l)M 


g(M,nx)« 


Y(1  +  nM)  -  (M2  -  1)^ 
1  +  (Y  -  1)M2 


For  a  skew  wire,  therefore,  we  can  re-define  the  sound  mode  sensitivity  coef¬ 
ficient  as: 

e^  ■  e^  +  A  cot  0  *  f(M)e^  +  [g(M,nx)  +  h(M,n^)  cot  (B-10) 

so  that  the  hot-wire  response  equation  becomes: 


-  e  a'  +  e_T'  +  e_  cot  0  s'  +  e  n ' 

i^T  T  T  s  rr 


(B-ll) 


Using  the  notation  of  Appendix  B,  Reference  5  we  can  write: 


e  cot  0 

T 


p  cot  0 
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so  that  Eqn.  (B-ll)  becomes: 


pr '  +  q?  '  +  ra '  +  srr 1 


CB- 12) 


or,  taking  the  mean  square  of  each  of  the  above  equation: 


2  , 2  ^  2  , 2  2  ,2  2  ,2 

a  T  +qf  +rcr  +•  s  n 


+  2pqt '5 '  +  2prT  'a  1 


(B-13) 


+  Ipst'tt'  +  2qrf  'a '  +  2qs£  n"  +  2rscj'n' 


Although  Eqn.  (B-13)  is  comprised  of  the  ten  unknowns: 


with  the  X-probe  we  can  only  find  6  unknowns  and  must  eliminate  4  unknowns  before 
proceeding.  As  in  Reference  1  we  rely  on  Morkovin's  suggestion  that  in  a  "good" 
wind  tunnel  the  sound  mode  is  uncorrelated  with  the  other  modes.  That  is,  we  let: 

T  V  “  ?  'tt"*'  ■  o'tt'  ™  o 
Furthermore  we  assume  that: 

■i 

c't'  ■  £(a ')2 • (t ' )2J  ,  i.e.,  that  R^T  ■  -  1.0 

With  the  above  assumptions,  Eqn.  (B-13)  reduces  to: 


2  ,2  ,  2  „  , ,  2  ,2  L  2,2 

p  T '  +q  *  +  r  a  +sn  + 

-  2pr(o,2)‘5(T,2)'-i  +  2qr?V 


2pqr'f 


(B-14) 


Letting  finally  p/r  ~  eT/e  s  x  and  dividing  both  sides  of  Eqn.  (B-14)  by  r' 
we  obtain:  _  c 

Y2  s  -  X2t'2  +  X2  cot20  5' 2  +  a'2 


4-  f(M)  +  [g(M,n  )  +  h(M,n  )  cot  0]  xl2  rr  * 4 

-  x  y  j 

+  2K2  cot  0  TV  -  2X  or  +  2X  cot  0  ?  'a' 


(B-15) 


which,  after  collecting  like  powers  of  X,  can  be  written  as: 

Y2  -  X2  £t'2  +  cot2  0  g,24*2cot  0  t'?'  +  (g  +  h  cot  0)2  rr'2J 

-  2X  j"crT  -  cot  0  5  'o 1  -  f(g  +hcot0)n'2 
r^2  .  c2  “721 


<B-  16) 


Y2  -  AX2  +  BX  +  C 


(B-17) 


For  an  X-probe  we  can  write  the  above  expression  for  both  wires  (denoted  1  and  2) 
and  modal  analysis  will  yield  the  following  quantities: 


T 1 2  +  cot2  0^  5'2  +2  cot  0^  T'5 '  +  (g  +  h  cot  0^)‘  tt'*  (B-18a) 


2  ,2 


1  -  -aT  +  cot  0^  5  'a 1  +  f<8  +  h  cot  0^)  n' 


1 2  ,  -2  f  2 

a  4*  f  tt 


<B-18b) 


(B* 18c) 


A2  -  T'2  +  cot2  02  g ,2  +  2  cot  02  T'C '  +  (  g  +  h  cot  02)2  n'2  (B-18d) 


B.  -  -aT  +  cot  09  ? 'aT  +  f  (g  +  h  cot  0„)  tt' 


(  B- 1 8  e  ) 


,2  r2  1 2 

«  +•  f  tt 


(B-18f) 


Note  that  must  equal  C2  and  therefore  Eqns.  (B-18c)  and  (B-18f)  reduce  to  a 
single  expression: 


■'2  +  f2  n'2 


( B- 18g) 


We  can  also  write  for  wires  1  and  2 


el'  “  P1T'  +  P1  cot  +  1'lCT'  +  ^Erl  +  (s  +  h  cot  p^]  n1  (B-19) 

e2'  -  P2t'  +  p2  cot  02  5  '  +  r2a'  +  [fr2  +  (g  +  h  cot  02>  p2J  rt 1  (B-20) 


Then,  taking  the  mean  square  of  the  product  'e  '  yields: 
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ra i ' «2  '  *  p1p2t'  +  plp2  COt  ?  It'*  “  P i r 2  aT 

+  p1p2  cot  0L  5ttT  +  P1P2  cot  0L  cot  02  $'2 

+  p^  cot  <t>x  J~'oV  -  P2rL  oT  +  P2r1  cot  02  ?  'o '  (B-21) 

+  rlr2  0  +  [ft^  +  (g  +  h  cot  0j)  pt]  [fr2 

+  (g  +  h  cot  02)  P2]  n'2 


e 


1 


Dividing  both  aides 


P1P2T'2  +  P1P2  COt  ^2  ^Tt'  +  P1P2  COt  ?T^' 

+  plp2  cot  0^  cot  02  ;'2  -  p^  at  +  PjT  cot  0L  FV 

’  P1P2  cT  +  P2rl  COt  ^2  * >a  *  +  rLr2  °'2  +  f2  flr2  11,2  ^B"22^ 
+  P^r2  f^8  +  h  c0‘‘  '&])  +  p2rl  f<8  +  h  cot  *52^rr'2 

+  P^P,p  (g  +  h  cot  0^)(g  *•  h  cot  0r) )  tt  ' 2 

of  (B-22)  by  r^r2  we  have! 


c1"  7 
1  2 

rir2 


Xi  [X2^T'2  +  cot  ? It'  +  cot  ? 'tI 

+  cot  0^  cot  02  ?'2  +  (g  +  h  cot  0^)(g  +  h  cot  02)  tt'2 

+  C-oT  +  cot  0^  ?'cr'  +  f  (g  +  h  cot  0^)  rr ' 2  j 

+  X2[-oT  +  cot  02  5 'a'  +  f(g  +  h  cot  0,;)  tt  * 2  ] 

+  [a'2  +  f2  n,23 

Xi  [D1X2  +  °2]  +  °3X2  +  °4 


3 

(B-23 ) 


(U-24) 


Equation  (B-24)  can  be  solved  by  the  "correlation  mode"  scheme  described  in 
Reference  5  to  yield: 
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(B-25a) 


D^  *  T*  +  cot  0^  cot  0^  5 '  +  (cor  0^  +  cot  02^  5 ' T 

+  (g  +  h  cot  0^ ) (g  +  h  cot  02)  tt'2 

D2  -  -ctT  +  cot  ;"'o  '  +  f  <g  +  h  cot  0L)  tt'2  ( B- 2 5 b ) 

D3  -  -ct  +  cot  02  5 ‘c '  +  f  <g  +  h  cot  02)  tt'2  (B-25c) 

D4  o'2  +  f2  n'2  (B*25d) 

Note  again  that  D,  »  C  and  that  Eqns.  (B»18b)  and  (B-18c)  are  identical  respec¬ 
tively,  to  Eqns.  (B-25b)  and  B-25c). 


Defining  A*  «  Bj 
to  the  following 

■  D2  and  tq  b2  ■  D3 1  the  governing  equations, 

expressions: 

therefore, 

,  reduce 

A1 

m 

T'2  +  cot20^ 

- 5  — r-r 

5 '  +  2cot  01  5 ' T k  +  (g  +  h 

cct  01)2 

(B-26a) 

A2 

m 

9  2 

T'  +  cot  02 

5 ,2+  2cot  02  5'IT~'  +  (g  +  h 

fN 

CM 

•Si 

4J 

O 

U 

~2 

(B-26b) 

A3 

m 

-<JT  +  cot  0^ 

?  'aT  +  f(g  +  h  cot  02)  n'2 

(B-26c) 

A4 

N 

-OT  +  COt  02 

+  f(g  +  h  cot  02>  tt'2 

(B-26d) 

D1 

m 

t'2  +  cot 

cot  02  5'2  +  (cot  0^  +  cot 

02)  fv 

(B-26e) 

+  (g  +  h  cot 

0j)  (s  +  h  cot  02)  tt ' 2 

C 

m 

~  M  c2  ~ 
a  +  f  tt' 

(  B-  2  6  f  ) 

Eqns.  (B-26a)  to  (B-26f)  represent  a  set  of  six  non-linear  simultaneous  equations 
which  can  be  solved  to  find  the  six  unknowns: 


and 


-109- 


APPENDIX  C 


CONE  DRAG  BUDGET 


(a)  Pressure  (Invlscld)  Drag 

The  cone  pressure  drag^p  is  made  up  of  the  curved-surface  drag  $c  and  the  base 
pressure  drag^t,,  The  former  is  the  product  of  the  axial  component  px  ■  p.  sin 

6  of  the  surface  pressure  pc,  times 


]  t 

i.  d 


the  cone  curved  area 

A  -  Jii- 
c  ain  0 


(C-l) 


Thus, 


A  “  pc  sln  0_li ^Te  "  pcnr 


(C-2) 


Furthermore  since  the  cone  is  slender  (  0  *  5°)  we  can  approximate  (for  -  3): 

2 

(03) 


P«u» 

2 

YM_ 


-la_ 

YM  2 


where  q  «  (1/2)  p  u  ,  the  dynamic  pressure.  Thus, 


#  -  -2a-  „r2 

C  YM„2 


The  base  pressure  can  only  be  estimated  from  the  available  data  sb  i 
b 


p  ^  I  p  .1  . 


4  "  4  YM.2 


and 


■0p  -  “Pb 


so  that 


(b)  Viscous  Pratt  (Laminar) 


Sines  the  laminar  ahaar  stress  tw  is  a  function  of  ths  distance  i  along  a  ray,  we 

have  to  integrate  its  x-component  twx 
■  tw  cos  0  over  the  cone  area: 

viscous  drag  <8  •» 

,  V  ,  (C-5) 

J  k  TwxdA  "  c°8  9  / 


From  the  sketch,  we  find  that  the 
surface  area  element  dA  is 


*  "  fsrf ttn  0  <c‘6> 


For  the  skin  stress  t  we  writs 

w 


2  PeUe  Cf  “  2  P«U*  Cf 


(C-7) 


where  cf  is  a  laminar  friction  coefficient: 

C  ^00800 

c,  ■  — — —  (where  Re  *  . -) 

f  /rT  8 


Since  s  ■  x/cos  0 


ct  - 


„  U 

- — -  /coa  9/x  (Re1  a  — ,  unit  Reynolds  No,) 


/rT1, 


(C-8) 


(C-9) 


Inserting  Eqna,  (C-6),  C-7),  and  (C-9)  into  (C-5)  we  got: 

AT'  '  cos  B  A 


^v  “  cos  0  J  J  “pl/coTA  “7I  “  dx  ckp 


tp  X 


_  Ztt  qc  tan  8  /cob  9  P 
/Re7  ** 

x“0 


4  qC  tan  9  f cos  6  ^3/2 

3n  AT 


-ill- 


Since  h  * 


t>  u«  1  Re 

- — r  and  Re'  ■  —  -  ■  — r~* 

2  tan  0  vw  t>  £> 

and  since  C  has  to  be  multiplied  by  /*3  to  apply  to  the  cone  (Reference  31,  p.1067) 
we  get 


r co3  a  •] 
o  J 


/Ret,  Vain  0 


(C-10) 


(c)  Viecc 


turbulent 


For  the  turbulent  case,  Eqn.  (C-8)  is  changed  to 


0.0262 


<RO 


^44 i 4 


2  » 


_ IS _ r  f5/7 

co177 


from  Reference  31,  p.  1109.  Here 


0.0262 


Re„  * 


us  u, 

00  to  X 


cos  0 


Since  0  is  small,  we  write  (cos  0)^7  !*  1.  Then 


kf5/7  -1/7 

Rf  ~<M7)  x  1 


and  T  ■  qc,  ta  Px 
w  f 


where  P  a 


(RO177 


Then  Eqn.  (C-5)  becomes 


Turbulent  drag  a  .gj 


and  since 


i’  r  -1/7 

J  J  P  tan  6  x  xd  xd0 

0  x 

hu/7 


h13/7p  -  *0—  t1« 

Re  U  '  ^4  tan  0J 
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<c-ll) 


with  Re  a  u  /v  ,  we  get 
V  00 


U  »  tan  E,)1" 

26  tan  0  Re, 


(d)  Definition  of  Drag  Coefficients 

If  the  total  cone  drag  35  it  defined  on  the  basia  of  the  base  area 


\ 

ao  that 


nr 


2 


$  "  3  +  •&  ■  qC  A  «  qA.(C  +  cb  > 

p  V  D  D  Pp 

we  can  compute  the  Cj,  for  each  contribution  as  followai 
Pressure  drag; 

From  Gqn.  (04)  we  get 


For  ■  3,  C  -  0.119. 

cn  p 

p 


Laminar  viscous  drag: 
From  Gqn .  (C-10)  we  get 


c„  ■ 

l  /Rep  / sin  8 

A  computation  will  be  made  for  the  lowest  pressure  experiment  (F  experiment, 
PQ  *  280  mm  Hg  abs.,  Rep  ■  69,000).  The  constant  C  is  taken  as  the  Blaaiuu 
0.664  reduced  by  107.  to  account  for  the  friction  decrease  due  to  Mach  3 
(over  the  incompressible  value).  We  find 

Cp  -  0.0248  (C-12) 


The  total  Cp  for  this  case  is 

Cft-C  +C  -0.144. 
5  \  ‘i 
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Turbulent  viscous  drag; 
From  Gqn.  (C-ll)  we  get 


C 

V 


T 


14  kf5/7  2  tan  5" 1 1 1 

13  tan  8  ..  ReD  J 


Here  the  computation  Is  done  for  the  A  experiment,  Re  ^  ■  178,000.  The 
reeult  Is 


C.  -  0.0284  (C- 13) 

dt 

The  total  C ^ for  this  case  Is 

ct>  “  cft  +(6  "  0.147  (C-14) 

^  J  m 


(e)  Comnarlaon  of  Invlscld,  Laminar  and  Turbulent  Drag 

We  note  from  Eqna .  (012)  and  013)  that  there  is  little  difference  between 
laminar  and  turbulent  friction  drag.  The  actual  difference  is  of  course  even 
smaller  because  In  the  A  configuration  only  a  portion  of  the  model  boundary- layer 
Is  turbulent.  This  near-equality  is  not  surprising,  even  accounting  for  the 
usual  higher  value  of  the  turbulent  friction  coefficient  over  the  laminar  one 
(Reference  31).  The  point  is  that  the  plots,  including  those  in  the  latter  Ref¬ 
erence,  show  the  flat  plate  coefficients  and  usually  at  higher  Rex  (from  105  to 
10').  In  the  present  case  the  Reynolds  numbers  are  lower,  and  the  laminar  fric¬ 
tion  coefficient  is  multiplied  by  the  factor  /3  to  account  for  the  conical  flow. 
Thus  the  present  laminar  coefficient  is  equal  to  or  larger  than  the  turbulent 
one  over  the  front  end  of  the  cone. 

Without  yet  looking  at  the  WEM-III  computer  results  which  give  the  "drag  dia¬ 
meter"  CqA,  we  ask  how  the  Cj%A  values  computed  above  compare  with  the  baBe  dia¬ 
meter  D.  For  the  A  test  condition  (the  post-BLT  case)  we  have 

C_A  -  (C  +  C  )A  -  0.147  “r-  -  0.115D2 

^  ^  1.  ^ 

(CdA)1|/2  -  0.34D,  0  -  2.94  (CQA)1/2 


with  negligible  difference  for  the  laminar  F  condition  (pre-BLt,  Eqn.  (C-12)). 
Thus,  6  base  diameters  (the  limit  of  the  present  test  surveys)  mean  a  distance 
of  17.7  (CqA)^ ,  However,  the  proper  similarity  parameter  is  the  wake  drag  dia¬ 
meter  (CpA)^,  Since  this  corresponds  only  to  the  viscous-friction  portion  C  A 
of  the  body  drag,  T 
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0.15D, 


(C^AjJ5  -  (C^A)1*  -  (0.0284  *£)*  - 

D  -  6.66  <CdA)w* 

for  the  turbulent  case  and,  similarly, 

-  °‘14D»  D  -  7.15  (CdA)w* 

for  the  laminar  case.  Thus  a  surveyed  distance  of  6D  equals  a  distance  of  about 
40  wake  drag  diameters  in  the  A  test  and  of  about  43  wake  drag  diameters  in  the 
F  test. 
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APPENDIX  D 


COMPUTER  PROGRAM  WEB- IV  (DRB4E) 


COMPUTATION  INSTRUCT IONS* 


7.1  Segregate  the  (Z,V)  from  the  (2,F)  points .  Note  that  the  first  point  of 

each  set  (i.e.,  each  current  code  A,  B,  etc.)  la  a  Z,  V  point.  The  points 
are  in  increasing  order  of  y  and  alternate  from  Z,  V  to  Z,  F. 


7.2  Assign  y's  to  each  Z,  V  or  Z,  F  point.  To  do  that)  first  note  the  first 

(y0)  and  laBt  (ye)  points  for  that  x;  see  Table  I.  All  current  codes  of  that 
x  obey  the  same  rules.  Then  add  0.010  inch  to  the  next  point  in  line,  all 
the  way,  until  the  last  point  ye. 


For  example,  the  card  on  p,  2  is  the  first  card  of  the  K  set,  x  ■  0.8,  WEMA, 


According  to  Table 
as  follows: 

l,  y0  ■  o.20o, 

ye  ■  0.595.  ThuB 

the  points  are 

listed 

Z.  V  Points 

Z 

.  F  Points 

.  y  2 

V 

_y 

Z 

F 

0.200  -.4722 

.2409 

0.205 

-.4769 

.4933 

0.210  -.4812 

.2406 

0.215 

-.4881 

.5132 

0.220  -.4967 

.2403 

0.225 

(first  point 

on  next 

•  • 

• 

* 

• 

• 

•  • 

s  • 

♦ 

• 

• 

s 

« 

* 

The  y  points  correspond  to  thoBe  of  the  WEM-I,  II  outputs  (where  they  are 
listed  as  y/D).  From  now  on,  we  work  with  the  Z,  V  points  separately  from 
the  Z,  F  points. 


7 . 3  Computation  of  the  "Sensitivity  Matrix"  (Z.  V  Data  Only) 

7.3.1  For  each  x,  go  to  Table  1,  select  the  "current  group".  Then  Table 
2  gives  the  current  1,  for  each  group,  matched  to  the  current  code, 
B,  C.....K.  At  each  y  of  that  particular  x, 


*  Numbering  of  subsections  in  this  Appendix  is  kept  identical  to  that  of  the 
actual  program  on  file  at  Philco-Ford.  Consistently,  references  to  page  and 
Table  numbers  refer  to  that  file,  not  to  page  or  Table  numbers  in  this  report. 


7.3.2 

7.3.3 
7.3.3a 

7.3.4 

7.3.5 

7.3.6 

7.3.7 

7.3.8 

7.3.9 

7.3.10 


Match  the  V  from  each  current  code  with  the  i  from  Table  2. 

You  now  have  an  i,  V  curve. 

Compute  R'  *  1000  V/i  for  each  1.  You  now  have  an  i ,  R'  curve. 
Compute  R  -  R'  -  1.42. 

Compute  by  least  squares  the  coefficients  C,  C,  and  C„  from  the 
curve  R,  i2!  it 


R  -  C  +  i2  +  C2  i4 

(That  is,  "plot"  R  versus  i2 
degree  polynomial  in  i2). 

Compute i 

,  i  1  i  dR  if 

w  2  R  di  R 

Compute: 


and  curve- fit  that  with  a  second' 


(C1  +  2C2  i2)  each  i 


Tw  "  <f~  ‘  l>  +  273.18  each  i 

o  o 

using  R  from  7.3.3  and  VQ,  RQ  from  Table  1.  (Note  that  which 
Ro  you  choose  depends  on  x) . 

Compute: 


aw 


o  o 
with  C  from  7.3.4, 


(f“  -  1)  +  273.18 


Compute : 
T 

wr 


T  -  T 
w  aw 

T 

aw 


from  3.3,7  and  3.3,6. 
Compute : 


K 


a  r 


w 


oo  R 


from  7.3.3  and  7.3.6  and  ®  ,  R  from  Table 

o  o 


Compute : 


each  i 


each  i 

1. 


c  .  -  0-992 _ 

cf  1  +  0.016  A  ' 

w 

from  7.3.5. 


each  i 
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7.3.11 

7. 

7. 

7. 

7. 


7. 

7.3.12 

7.3.13 


Obtain  (at  the  y  considered)  the  quantity  Re  ■  2  (Wire  Reynolds 

No.) 


>  .11a 


l.llb 


At  each  point  y  compute: 

k  -  (6.47  +  0.02268  (T 


aw 


390))  x  10 


-4 


with  T  from  step  7.3.7, 
aw 


Compute: 

Nu 

in 


o^o  1  C  .  -• 6 

~  t  57  x  10 


where  ot  Rq  are  from  Table  1,  n  *  3.14......  k  is  from 

7.3.11a,  C  and  from  7.3.4  and  K  (in  inches)  is  also  given 
from  Table  1. 


.11c  At  each  y,  we  now  have  a  value  of  Num  from  7.3.11b  and  a 
value  of  Re  as  given  in  7.3.11.  Using  the  points  Num,  Re 
for  the  particular  x-Station  only,  curve-fit  these  points 
into  a  single  straight  line: 

Nu  ■  N.  +  N.,  Re 
m  12 

and  thus  find  N^  and  N^ • 

.lid  Similarly,  curve-fit  the  ratio  Taw/T0  versus  Re: 

T  ? 

T]  is  »  Hl  +  H2(Re)  +  H^(Re)^ 

0 

where  Taw  is  from  step  7.3.7  and  T  is  from  the  WEM-1I  out¬ 
put,  designated  TOTAL  TEMP  (deg.  K) .  Thus  find  the  constants 
Hi,  H2  and  H3. 

.He  Print  the  quantities  Ni,  N2,  Hi,  H2  and  H3  on  top  of  the  out¬ 
put  page  called  "INTENSITIES"  for  that  particular  x-Station. 


Compute : 

0  "  Nu  N2  “  N,  +  N  Re 

m  12 

with  Re  from  7 .3.11. 


Compute : 

*  -  {p  <H2  +  2H3  Re) 

from  Re  from  7.3.11,  T|  *  T  /T  ,  H„  and  H,.  from  7. 3. lid. 

1  aw  o  2  J 
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7.3.14  Compute! 


: .3.15 


e  *  C-  A  - sL 

m  f  w  l.  t  _ 

wr 


each  1 


using  C £  from  7.3.10,  A^'  from  7.3.5,  r>Air  from  7.3.8,  and  0,  \|i 
from  7.3.12  and  7.3.13. 


Compute : 


r  a  ' 

Cf  K  +  A  '  (K  +  0.7650  -  1.885)  -  0.765\l^“~J 


each  i 


with  Cf,  Ay',  0,  ill ,  and  Twr  as  before  and  K  from  7.3.9, 


,.3.16  Compute: 


0  0.41  M 


1  +  0.2  M 


where  M  is  the  Mach  number  at  this  particular  y  (from  the  output 
of  WEM- II). 


7.3.17  Compute: 


Be„,  -  e 
T  m 

le_  +  e 
T  m 


each  i 


using  7.3.14,  7.3.15  and  7.3.16.  Note :  It  is  desired  to  print 
out  the  matrix  i,  y,  X(the  "sensitivity  matrix"). 

Experiment:  .  X  *  _ 

_ i  corresponding  to  current  codes  (see  Table  2) _ _ 

y  Code  A  Code  3  Code  C  Code  D  . K. 


.  Values  of  X  from  7.3.17. 

7.3.1H  Print  the  "Wire  Voltage  Matrix"  for  the  V  data: 

1  corresponding  to  current  codes  (see  Table  2) 
Y/D  Code  A  Code  B  ... 


.  Values  of  V  at  each  point,  fron  input  cards. 

Print  this  between  the  "sensitivity  matrix"  and  the  "mode  matrix" 
for  the  Z,  V  data. 
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7 .4  Computation  of  the  Mode  Matrix  for  Wideband  Fluctuations 


At  each  x,  consider  now  the  points  which  have  Z,  V  data  only.  At  eacy  y; 
7.4.1  Match  the  Z  and  V  from  each  current  code  with  the  i  from  Table  2. 
7. 4. la  Compute: 

Z]_  "  each  i 


7.4.2 


7,4.3 


where  ZN,  tj_  come  from  Table  1  (one  for  each  experiment  and 
x-Station) . 


2  2 

If  (Z  -  )  is  negative,  set  to  zero. 


Compute : 

Z' 

where  G  comes  from  Table  1. 


21 

GV 


each  i 


Compute : 

..2 


‘  2'  T 

-*eT  +  V 


each  i 


using  Z'  from  7.4.2,  a  from  7.3.16,  eT  from  7.3.15  and  e  from 
7.3.14.  1  m 

2 

7.4.4  For  each  i,  we  now  have  an  X(from  7.3.17)  and  a  Y  (from  7.4.3). 
Solve,  by  curve-fitting,  a  second  degree  polynomial  (least  squares) 
for  the  coefficients  A,  B,  C  of 

Y2  -  AX2  +  BX  +  C  each  y 

(do  not  confuse  this  with  the  C  of  step  7.3.4)! 

7.4.5  Call: 

T  -  (A) 4 

a  m  (Cf11  each  y 

r  _  *  B/2  <tt 

<TT 


7.4.6 


so  that  at  each  y  we  have  a  t,  a  a,  and  an  r  , 

aT 

Compute : 

SRA  *  - —  each  y 

0.4  T 
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Note!  Print  out  the  X-Y  matrix! 


Consider  now  Che  y  point  (at  each  x)  having  Z,  F  data  only.  For  each  y, 


HP!’. 1  ■ 


7.5.7 


7.5.8 


Again  form: 
Y  ^  M 

TF 


<“'t  *  %>2  v2 


each  i 


as  in  atep  7.5.3. 

Repeat  the  procedure  of  step  7.5.4:  Find  App,  Bpp,  Cpp  from: 


FF 


A__X  B__X  +  C__ 

FF  FF  FF 


with  X  from  7.3.17  as  always. 
7.5.9  Form: 

A  m  al  tm. 

u  4  aff 


7.5.10 


with  u'  ■  63(600  u  (us  is  from  WEM-II  output).  k„  is  from  step 
7.5.4,  App  from  step  7.5.8. 

Also  form: 

A  u'  C?_ 

At  "  r  cpp 


7.5.11 


with  u',  Cp,  CpF  as  per  above. 
Also  form: 


Au  Au 
D  "  2.54 

where  Au  is  from  step  7.5.9  (i.e.,  column  LAMBDA (U)  of  printout) 
and 


A 


_J _ 

2.54 


with  A^  from  step  7.5.10  (i.e.,  column  LAMBDA(T)  of  printout). 


Note: 


Print  the  matrix  for  the  X,  i,  Yp,  and  X,  i  and  Y_p  in  a  way 
exactly  analogous  to  the  matrix  *  X,  i,  Y  of  pagetl  12. 


rr: 
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-a.aai  J  a.\,.aa-i .  i A  1 


-1l<l  l-l  lltwi'.'tuilt  :lf|i  l-n  .• 


FINAL  PRINTOUT 


INTENSITIES 


Experiment i 


Y/D  TAU  SIGMA  R(ST) 


■u-,7.u.g 


u  1  /Au  T7AT 


SCALES 


Y/D  LAMBDA(IJ)  LAMB  DA  (T 'l 


LAMOn/D  LAM(TWD 


** i*w**V«  '•  tt-iiw,*  J.  1  "  ■’* 


APPENDIX  E 


A  CRITERION  ESTABLISHING  A  MINIMUN  REYNOLDS  NUMBER  FOR  THE 
EXISTENCE  OF  TURBULENCE 


It  would  be  helpful  to  devise  a  rule  establishing  a  state  of  fluid  flow  below 
which  turbulence  cannot  "exist".  There  is  presently  controversy  as  to  whether 
the  turbulent  boundary- layer  remnant  shed  into  a  near  wake  continues  "seeding" 
the  near  wake  with  turbulence,  or  whether  it  get  completely  and  quickly  obliter¬ 
ated  by  the  shoulder  expansion.  Theoretical  estimates  of  near-wake  turbulence 
behavior  (References  18  and  32)  already  imply  that  the  effectiveness  of  such  anti- 
transition  hinges  oh  the  magnitude  of  the  turbulence  Reynolds  number 

Re.  - 

A  v 

where  u’  is  the  velocity  fluctuation,  A  the  turbulence  scale  and  v  the  kinematic 
viscosity . 

A  very  approximate  criterion  can  be  found  as  follows:  a  fluid  element  moving 
with  speed  u'  relative  to  its  adjacent  fluid  loses  this  speed  because  of  viscous 
friction.  On  purely  theoretical  grounds  the  element  has  a  finite  lifetime  t* 
beyond  which  its  speed  finally  becomes  uniform  with  that  of  the  adjacent  fluid, 
and  thus  the  element  ceases  being  identifiable.  During  its  lifetime  Tt  each 
element  traverses  a  "life-distance"  s  ■  u'/%.  To  qualify  as  a  component  of  tur¬ 
bulent  flow,  an  element  must  have  a  finite  life-time  Tt  and  a  distinguishable 
life-distance  s.  Here  we  postulate  that  the  minimum  s  for  element  identification 
is  the  element  ("eddy")  size  itself  (i.e.,  turbulence  1b  not  defined  to  exist 
unless  the  eddies  travel  a  distance  at  least  equal  to  their  own  size);  that  Is, 
for  turbulence  to  exist  we  require 

u't  ■  e  ■  A  (E- 1 ) 

Kovasznay  (Reference  4)  finds  for  the  lifetime  of  eonvected  vorticlty: 


k  v  (2tt)‘ 


(E-2) 


so  that  when  we  combine  (E-l)  and  (E-2)  we  get  a  turbulence  Reynolds  number 


>  MiA  «  . 


(2rr)  -  40 


below  which  "turbulence  cannot  exist". 
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W  IPONIOHINQ  MIU  T  ANi  A  C  T  I  V  I  T  V 


i  A  m  I  »«  A  C  T 


The  study  of  turbulence  on  a  flat  plate  at  Mach  9.4  has  been  completed  with  the  dis¬ 
covery  that  the  pressure  fluctuations  predominate,  that  the  velocity  fluctuations  are 
very  small  and  that  wall-cooling  has  a  strong  influence  in  decreasing  the  temperature 
fluctuations.  Tests  have  been  performed  at  AEDC/B  to  extend  these  results  to  the  cone 
geometry,  to  broaden  the  Mach-  and  Reynolds-numbevs  matrix  and  to  investigate  ablation 
effects.  Although  still  incomplete,  analysis  of  the  data  Bhow  very  large  temperature 
fluctuationa  due  to  the  adiabatic  wall,  relatively  little  change  of  the  turbulence  level 
past  the  transition  point,  and  discernible  but  small  fluctuation  increase  for  small 
mass  injection  rates.  Correlations  have  began  of  all  the  available  boundary-layer  tur 
bulence  data  for  the  purpose  of  scaling  the  fluctuation  distribution  to  arbitrary  flight 
conditions.  Thus  far  the  available  correlation  schemes  fit  the  velocity  fluctuation 
data  only  in  the  outer  half  of  the  layer.  Better  success  is  enjoyed  for  the  temperature 
difference  between  the  wall  and  the  external  flow.  Experiments  were  also  conducted  to 
measure  the  traverse  velocity  fluctuation  component  and  the  turbulent  shear  streaaes 
with  specially-developed  hot-wire  probes;  these  data  are  now  available  on  tape.  In  the 
near-wake  investigation  the  completed  mean-flow  measurements  were  analyzed  and  Bhow  that 
a  "plateau"  tends  to  form  as  the  Reynolds  number  increases  beyond  BLT.  The  turbulent 
viscosity  was  measured  and  found  to  decrease  beyond  BLT  as  predicted  by  the  AVCO  theor 
Other  theoretical  predictions,  Buch  as  the  radial  distribution  of  velocity  beyond  the 
neck,  were  also  verified.  Fluctuation  measurements  revealed  that  the  shoulder  expan¬ 
sion  nearly  destroys  the  velocity  but  not  the  temperature  fluctuations,  and  that  the 
wake  is  "double-structured"  beyond  the  neck,  The  study  of  electron,  velocity  and  tem 
oarature  soaca-time  correlations  in  the  far  WAke  has  been  completed. 
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